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Abstract 



A perfect number is a positive integer N such that the sum of all the pos- 
itive divisors of N equals 2N, denoted by cr(N) = 2N. The question of the 
existence of odd perfect numbers (OPNs) is one of the longest unsolved problems 
of number theory. This thesis presents some of the old as well as new approaches 
to solving the OPN Problem. In particular, a conjecture predicting an injec- 
tive and surjective mapping X = a ^ k - , Y = - between OPNs N = p k m 2 
(with Euler factor p k ) and rational points on the hyperbolic arc XY = 2 with 
1<X<1.25<1.6<F<2 and 2.85 < X + Y < 3, is disproved. Various results 
on the abundancy index and solitary numbers are used in the disproof. Numeri- 
cal evidence against the said conjecture will likewise be discussed. We will show 
that if an OPN N has the form above, then p k < |m 2 follows from [15]. We will 
also attempt to prove a conjectured improvement of this last result to p h < m by 
observing that i an d 21eJ. I^L i n all cases. Lastly, we also prove the fol- 

r 

lowing generalization: If N = Y\ Vi ai is the canonical factorization of an OPN N, 

i=l 

then a(pi ai ) < for all i. This gives rise to the inequality N 2 ~ r < (|)(f) r_1 , 
which is true for all r, where r = u(N) is the number of distinct prime factors of N. 



Chapter 1 

The Problem and Its Background 



1.1 Introduction 

Number theory is that branch of pure mathematics concerned with the properties 
of integers. It consists of various results and open problems that are easily under- 
stood, even by non-mathematicians. More generally, the field has developed to 
a point where it could now tackle wider classes of problems that arise naturally 
from the study of integers. 

A particular example of an unsolved problem in number theory, one which 
has long captured the interest of both amateur and professional mathematicians, 
is to determine whether an odd perfect number exists. A positive integer n which 
is equal to the sum of its positive proper divisors (that is, excluding n itself) is 
called a perfect number. The smallest known example of a perfect number is 6 
since 1 + 2 + 3 = 6, where 1, 2 and 3 are the positive proper divisors of 6. The 
Pythagoreans considered such numbers to possess mystical properties, thus calling 
them perfect numbers. 



The sum of all positive divisors of a positive integer n is called the sigma 
function of n, denoted by a(n). The definition of a perfect number is thus equiv- 
alent to determining those n for which u{n) —n = n. We formalize this definition 
as follows: 

Definition 1.1.1. A positive integer TV is perfect if cr(N) = 2N. 
For example, cr(6) = 1 + 2 + 3 + 6 = 12 = 2(6). 

Even though the Greeks knew of only four perfect numbers, it was Euclid 
who proved that if the sum 1 + 2 + 2 2 + . . . + 2 P ~ 2 + 2 P ~ 1 = TP — 1 is prime, then 
2 P_1 (2 P — 1) is perfect. Consider the sum 1 + 2 = 3 as an example. Since 3 is 
prime, then 2 2_1 (2 2 — 1) = 2(3) = 6 is perfect. 

Euler subsequently proved the following theorem about (even) perfect num- 
bers: 

Theorem 1.1.1. Every even perfect number is of the form N = 2 P ~ 1 (2 P — 1), 
where P and 2 P — 1 are primes. 

It is easy to prove that if C is composite, then 2 C — 1 is also composite. 
However, if P is prime, it does not necessarily follow that 2 P — 1 is also prime. 
(Consider the case P = 11 : 2 11 — 1 = 2047 = 23 • 89, which is composite.) 
Primes of the form 2 P — 1 are called Mersenne primes, after the French monk 
and mathematician Marin Mersenne. In view of Theorem 11.1. 11 the problem of 
searching for even perfect numbers is thus reduced to looking for Mersenne primes, 
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since the theorem essentially says that the even perfect numbers are in one-to-one 
correspondence with the Mersenne primes. 

As of this writing, 44 perfect numbers are known [70], the last few of which 
have been found with the aid of high-speed computers. Following is a list of all 
the known exponents p for which M p = 2 P — 1 is prime, along with other pertinent 
information (note that P p refers to the perfect number N = 2 P ~ 1 (2 P — 1)): 



number 


p (exponent) 


digits in M p 


digits in P„ 


year 


discoverer 


1 


2 


1 


1 






2 


3 


1 


2 






3 


5 


2 


3 






4 


7 


3 


4 






5 


13 


4 


8 


1456 


anonymous 


6 


17 


6 


10 


1588 


Cataldi 


7 


19 


6 


12 


1588 


Cataldi 


8 


31 


10 


19 


1772 


Euler 


9 


61 


19 


37 


1883 


Pervushin 


10 


89 


27 


54 


1911 


Powers 


11 


107 


33 


65 


1914 


Powers 


12 


127 


39 


77 


1876 


Lucas 


13 


521 


157 


314 


1952 


Robinson 


14 


607 


183 


366 


1952 


Robinson 


15 


1279 


386 


770 


1952 


Robinson 



number 


p (exponent) 


digits in M p 


digits in P p 


year 


discoverer 


16 


2203 


664 


1327 


1952 


Robinson 


17 


2281 


687 


1373 


1952 


Robinson 


18 


3217 


969 


1937 


1957 


Riesel 


19 


4253 


1281 


2561 


1961 


Hurwitz 


20 


4423 


1332 


2663 


1961 


Hurwitz 


21 


9689 


2917 


5834 


1963 


Gillies 


22 


9941 


2993 


5985 


1963 


Gillies 


23 


11213 


3376 


6751 


1963 


Gillies 


24 


19937 


6002 


12003 


1971 


Tuckerman 


25 


21701 


6533 


13066 


1978 


Noll, Nickel 


26 


23209 


6987 


13973 


1979 


Noll 


27 


44497 


13395 


26790 


1979 


Nelson, Slowinski 


28 


86243 


25962 


51924 


1982 


Slowinski 


29 


110503 


33265 


66530 


1988 


Colquitt, Welsh 


30 


132049 


39751 


79502 


1983 


Slowinski 


31 


216091 


65050 


130100 


1985 


Slowinski 


32 


756839 


227832 


455663 


1992 


Slowinski, Gage et al. 


33 


859433 


258716 


517430 


1994 


Slowinski, Gage 


34 


1257787 


378632 


757263 


1996 


Slowinski, Gage 


35 


1398269 


420921 


841842 


1996 


Armengaud, Woltman, et al. 



number 


p (exponent) 


digits in M p 


digits in P p 


year 


discoverer 


36 


2976221 


895932 


1791864 


1997 


Spence, Woltman, et al. 


37 


3021377 


909526 


1819050 


1998 


Clarkson,Woltman, et al. 


38 


6972593 


2098960 


4197919 


1999 


Hajratwala, Woltman, et al. 


39 


13466917 


4053946 


8107892 


2001 


Cameron, Woltman, et al. 


?? 


20996011 


6320430 


12640858 


2003 


Shafer, Woltman, et al. 


?? 


24036583 


7235733 


14471465 


2004 


Findley, Woltman, et al. 


?? 


25964951 


7816230 


15632458 


2005 


Nowak, Woltman, et al. 


?? 


30402457 


9152052 


18304103 


2005 


Cooper, Woltman, et al. 


?? 


32582657 


9808358 


19616714 


2006 


Cooper, Woltman, et al. 



Question marks (??) were used instead of a number for the the last of the Mersenne 
primes because it will not be known if there are other Mersenne primes in between 
these until a check and double check is done by the Great Internet Mersenne Prime 
Search (GIMPS) [70] or other similar coordinated computing projects. 

Note that all of the known perfect numbers are even. 
1.1.1 Statement of the Problem 

Our primary object of interest would be odd perfect numbers, though some of 
the results that would be discussed apply to even perfect numbers as well. It is 
unknown whether there are any odd perfect numbers. Various results have been 
obtained, but none has helped to locate one or otherwise resolve the question of 
their existence. 
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This thesis explores some of the old as well as new approaches used in trying 
to solve the Odd Perfect Number (OPN) Problem, namely: 

• the use of the abundancy index to derive conditions for the existence of odd 
perfect numbers, in such instances as: 

— bounding the prime factors of an OPN; 

— determining whether a particular rational number may be an abun- 
dancy index of a positive integer; 

— increasing the lower bound for the number of distinct prime factors, 
u(N), that an OPN N must have; 

• proving the inequality p k < |m 2 where N = p k m 2 is an OPN with p k the 
Euler factor of N, p = k = 1 (mod 4), and gcd(j>, m) — 1; 

• attempting to prove the conjectured improvement of the result p k < |m 2 to 
p k < m by observing that ^ 1 and ^ apply in all cases; 

r 

• generalizing the result p k < |m 2 to: if N — Y\_Pi ai is the canonical factor- 

i=l 

ization of an OPN N, then a{pi a *) < for all i; 

• showing that N 2 ~ r < (|)(|) r ~ 1 follows from <r(pi ai ) < for all i, where 
r = u(N) is the number of distinct prime factors of an OPN N; 

• disproving the conjectured injectivity and surjectivity of the mapping 

X = 4p-, Y = between OPNs N = p k m 2 (with Euler factor p k ) and 
rational points on the hyperbolic arc XY = 2 with 
1 < X < 1.25 < 1.6 < Y < 2 < 2.85 < X + Y < 3; 



• a host of other interesting results on perfect numbers, including: 

— establishing that two consecutive positive integers cannot be both per- 
fect; 

— the proof that an OPN is a sum of two squares. 

Specifically, this thesis, among other things, aims to present (partial) expositions 
of the following articles/notes that deal with some of the above concerns: 

• u Abundancy 'Outlaws' of the Form a ^ +t n by W. G. Stanton, which was a 
joint research project with J. A. Holdener at Kenyon College, Gambier, OH 
(2007); 

• " Conditions Equivalent to the Existence of Odd Perfect Numbers" by J. A. 
Holdener, which appeared in Mathematics Magazine 79(5) (2006); 

• 11 Bounding the Prime Factors of Odd Perfect Numbers", an undergraduate 
paper of C. Greathouse at Miami University (2005); 

• " Hunting Odd Perfect Numbers: Quarks or Snarks?" by J. McCleary of Vas- 
sar College, Poughkeepsie, NY, which consists of lecture notes first presented 
as a seminar to students of Union College, Schenectady, NY (2001); 

• " Consecutive Perfect Numbers (actually, the Lack Thereof!)" by J. Riggs, 
which was an undergraduate research project with J. A. Holdener at Kenyon 
College (1998). 



1.1.2 Review of Related Literature 

Benjamin Peirce was the first to prove (in 1832) that an OPN N must have at 
least four distinct prime factors (denoted u(N) > 4) |53j. Seemingly unaware of 
this result, James Joseph Sylvester published a paper on the same result in 1888, 
hoping that research along these lines would pave the way for a general proof 
of the nonexistence of an OPN [6S]. Later that same year, Sylvester established 
u)(N) > 5 [69]. This marked the beginning of the modern era of research on 
OPNs. It was only in 1925 that Gradstein was able to improve Sylvester's result 
to u(N) > 6 [19J . In the early 1970s, Robbins and Pomerance independently 
established that u(N) > 7 [55J. Then, Chein demonstrated that u(N) > 8 in his 
1979 doctoral thesis [6], which was verified independently by Hagis in a 1980 200- 
page manuscript [23]. Recently, in a 2006 preprint titled "Odd perfect numbers 
have at least nine distinct prime factors" , Pace Nielsen was able to prove that 
u(N) > 9, and if N 2 is an OPN which is not divisible by 3, that co(N 2 ) > 12. "The 
proof ultimately avoids previous computational results for odd perfect numbers." 



Sylvester also showed (later in the year 1888) that an OPN cannot be divis- 
ible by 105 [69] , and Servais (in the same year) proved that the least prime divisor 
of a perfect number with r distinct prime factors is bounded above by r + 1 [62] . 
Grim [22], Cohen and Hendy [9], and McDaniel [49] established improvements 
and extensions to this last result later on. 



Dickson showed in 1913 that there can only be finitely many OPNs with r 
distinct prime factors, for a given positive integer r [16]. Kanold then generalized 
Dickson's theorem in 1956 to include any positive integer n satisfying = 2 ( 
where a and b are positive integers and b ^ 0. [43] 

(Call a number n non-deficient if > 2 . Dickson called a number prim- 
itive non- deficient provided that it is not a multiple of a smaller non-deficient 
number.) Mathematicians considered Dickson's approach (i.e. first delineating 
all of the finitely many primitive odd non-deficient numbers associated with a 
particular r-value and then determining which among them are equal to the sum 
of their positive proper divisors) to the OPN question to be impractical for most 
values of r, making it necessary to explore alternative approaches to examining 
the possible structure of an OPN. Pomerance suggested the following class of the- 
orems in 1974: An OPN is divisible by j distinct primes > N [55]. Kanold was 
successful with j — 1, N — 60 in 1949 and used only elementary techniques [41] . 
In 1973, with the aid of computation, Hagis and McDaniel improved Kanold's 
finding to j = 1, N = 11200 [25]. This was pushed to j = 1, N = 100110 by the 
same authors in 1975 [26]. Pomerance showed that j = 2, N = 138 in the same 
year. [56] 

In a recent preprint (titled 11 Odd perfect numbers have a prime factor ex- 
ceeding 10 8 ") that appeared 2006, authors Takeshi Goto and Yasuo Ohno report 
that the largest prime factor of an OPN exceeds 10 8 [IS] . It betters the previous 
bound of 10 7 established by Jenkins in 2003 [39] . New bounds for the second and 
third largest prime factors of an OPN were given by Iannucci in 1999 [36] and 2000 
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when he announced that they are larger than 10 and 10 2 , respectively. 



Mathematicians also began considering estimates on the overall magnitude 
of an OPN by imposing lower bounds. Turcaninov obtained the classical lower 
bound of 2 ■ 10 6 in 1908 [5j. The following table summarizes the development of 
ever-higher bounds for the smallest possible odd perfect number: 



Author 


Bound 


Kanold (1957) 


10 20 


Tuckerman (1973) 


10 36 


Hagis (1973) 


1Q 50 


Brent and Cohen (1989) 


10 160 


Brent et al. (1991) 


10 300 



There is a project underway at http:/ /www. oddperfect.org (organized by William 
Lipp) seeking to extend the bound beyond 10 300 . A proof for 10 500 is expected 
very soon, as all the remaining factorizations required to show this are considered 
"easy", by Lipp's standards. 



It would also be possible to derive upper bounds on the overall size of an 
OPN in terms of the number of its distinct prime factors. Heath-Brown was able 
to show, in 1994, that if n is an odd number with cr(n) = an, then n < (4d) 4 , 
where d is the denominator of a and r is the number of distinct prime factors 
of n |30j. Specifically, this means that for an OPN n, n < 4 4r which sharpens 

9 

Pomerance's previous estimate of n < (4ry in 1977 [56]. Referring to his 
own finding, Heath-Brown remarked that it still is too big to be of practical value. 



Nonetheless, it is to be noted that if it is viewed alongside the lower bound of 
10 300 given by Brent et. al. [3], then Sylvester's 1888 result that u(n) > 5 could 
then be demonstrated by no longer than a footnote. 3 In 1999, Cook enhanced 
Heath-Brown's result for an OPN with r distinct prime factors to n < D 4V where 
D = (195) 1/7 « 2.124 pi]. In 2003, Pace Nielsen refined Cook's bound to n < 2 4r 

Addressing the OPN question from a congruence perspective on the allow- 
able exponents for the non-Euler prime factors, Steuerwald showed in 1937 that 
if 

n = p * qi Wi q2 W2 . .. q Ws 

was an OPN where p, q±, q%, ■ ■ ■ , q s are distinct odd primes and p = a = 1 (mod 4), 
then not all of the /3j's can equal 1 [55] ■ Further, Kanold discovered in 1941 that 
it is neither possible for all /3j's to equal 2 nor for one of the /Vs to be equal to 2 
while all the rest are equal to 1 |4T)] . Hagis and McDaniel proved in 1972 that not 
all the /3j's can be equal to 3 [23]. Then in 1985, Cohen and Williams summarized 
all previous work done on this area by eliminating various possibilities for the /3j's, 
on the assumption that either some or all of the /3j's are the same [TU] . 

In 2003, Iannucci and Sorli placed restrictions on the /3j's in order to show 
that 3 cannot divide an OPN if, for all i, fli = 1 (mod 3) or = 2 (mod 5). They 
also provided a slightly different analysis by giving a lower bound of 37 on the 
total number of prime divisors (counting multiplicities) that an OPN must have 
HO 300 < n < 4 4r implies that r > 4.48. 



s s 

(i.e. they proved that if n = p a \\ q l Wl is an OPN, then Q(n) = a + 2^ ft > 37) 

i=l i=l 

|38j . This was extended by Hare later in the year 2003 to Q(n) > 47 [28J. In 
2005, Hare submitted the preprint titled "New techniques for bounds on the total 
number of prime factors of an odd perfect number" to the journal Mathematics of 
Computation for publication, where he announced a proof for Q(n) > 75 [25] . 



In order to successfully search for perfect numbers, it was found necessary 

to consider a rather interesting quantity called the abundancy index or abundancy 

ratio of n, defined to be the quotient ° W . Obviously, a number n is perfect if 

n 

and only if its abundancy index is 2. Numbers for which this ratio is greater than 
(less than) 2 are called abundant {deficient) numbers. 



It can be shown that the abundancy index takes on arbitrarily large values, 
we can make the abundancy index to be as close to 1 as we please because 

7) -I- 1 

for all primes p. In fact, Laatsch showed in 1986 that the set 



p p 

of abundancy indices °" U for n > 1 is dense in the interval (l,oo) H31. (Let 

n 

I(n) = a , and call a rational number greater than 1 an abundancy outlaw if it 
n 

fails to be in the image of the function I. [34J) Interestingly, Weiner proved that 
the set of abundancy outlaws is also dense in (l,oo)! [73] It appears then that 
the implicit scenarios for abundancy indices and outlaws are both complex and 
interesting. 



In 2006, Cruz [12J completed his M. S. thesis titled "Searching for Odd 
Perfect Numbers" which contained an exposition of the results of Heath-Brown 
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[3D] and Iannucci/Sorli [38J. Cruz also proposed a hypothesis that may lead to a 
disproof of the existence of OPNs. 



Chapter 2 

Preliminary Concepts 



The concept of divisibility plays a central role in that branch of pure mathematics 
called the theory of numbers. Indeed, mathematicians have used divisibility and 
the concept of unique factorization to establish deep algebraic results in number 
theory and related fields where it is applied. In this chapter, we survey some 
basic concepts from elementary number theory, and use these ideas to derive the 
possible forms for even and odd perfect numbers. 

2.1 Concepts from Elementary Number Theory 

For a better understanding of the topics presented in this thesis, we recall the 
following concepts. 

Definition 2.1.1. An integer n is said to be divisible by a nonzero integer m, 
denoted by m \ n, if there exists some integer k such that n = km. The notation 
m \ n is used to indicate that n is not divisible by m. 

For example, 143 is divisible by 11 since 143 = 11-13. In this case, we also 
say that 11 and 13 are divisors /factors of 143, and that 143 is a multiple of 11 
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(and of 13). On the other hand, 143 is not divisible by 3 since we will not be able 
to find an integer k that will make the equation 3k = 143 true. 

If n is divisible by m, then we also say that m divides n. 

We list down several properties of divisibility in Theorem 2.1.1. 
Theorem 2.1.1. For integers k, I, m, and n, the following are true: 

• n | ; 1 | n, and n \ n. (Any integer is a divisor of 0, 1 is a divisor of any 
integer, and any integer has itself as a divisor.) 

• m | 1 if and only if m = ±1. (The only divisors of 1 are itself and —1.) 

• If k\m and I \ n, then kl \ mn. (Note that this statement is one-sided.) 

• If k | / and I \ m, then k \ m. (This means that divisibility is transitive.) 

• m | n and n \ m if and only if m = ±n. (Two integers which divide each 
other can only differ by a factor of ±1.) 

• If m | n and n ^ ; then \m\ < \n\. (If the multiple of an integer is nonzero, 
then the multiple has bigger absolute value than the integer.) 

• If k | m and k \ n, then k \ (am + bn) for any integers a and b. (If an 
integer divides two other integers, then the first integer divides any linear 
combination of the second and the third.) 



A very useful concept in the theory of numbers is that of the GCF or GCD 
of two integers. 

Definition 2.1.2. Let m and n be any given integers such that at least one of 
them is not zero. The greatest common divisor of m and n, denoted by gcd(m, n), 
is the positive integer k which satisfies the following properties: 

• k | m and k \ n; and 

• If j \ m and j \ n, then j | k. 

Example 2.1.1. The positive divisors of 36 are 1, 2, 3, 4, 6, 9, 12, 18 and 36. 
For 81 they are 1, 3, 9, 27 and 81. Thus, the positive divisors common to 36 and 
81 are 1, 3 and 9. Since 9 is the largest among the common divisors of 36 and 81, 
then gcd(36,81) = 9. 

Another concept of great utility is that of two integers being relatively prime. 

Definition 2.1.3. Let m and n be any integers. If gcd(m,n) = 1, then m and n 
are said to be relatively prime, or coprime. 

Example 2.1.2. Any two consecutive integers (like 17 and 16, or 8 and 9) are 
relatively prime. Note that the two consecutive integers are of opposite parity 
(i.e. one is odd, the other is even). If two integers are of opposite parity, but 
not consecutive, it does not necessarily follow that they are relatively prime. (See 
Example [2XU) 



It turns out that the concept of divisibility can be used to partition the set 
of positive integers into three classes: the unit 1, primes and composites. 



Definition 2.1.4. An integer P > 1 is called a prime number, simply a prime, if 
it has no more positive divisors other than 1 and P. An integer greater than 1 is 
called a composite number, simply a composite, if it is not a prime. 

There are only 23 primes in the range from 1 to 100, as compared to 76 
composites in the same range. Some examples of primes in this range include 2, 
7, 23, 31, 41, and 47. The composites in the same range include all the larger 
multiples of the aforementioned primes, as well as product combinations of two 
or more primes from the range 1 to 10. (Note that we get 100 by multiplying the 
two composites 10 and 10.) We casually remark that the integer 2 is the only 
even prime. The integer 1, by definition, is neither prime nor composite. We shall 
casually call 1 the unit. 

From the preceding discussion, we see that the set of prime numbers is not 
closed with respect to multiplication, in the sense that multiplying two prime 
numbers gives you a composite. On the other hand, the set of composite numbers 
is closed under multiplication. On further thought, one can show that both sets are 
not closed under addition. (It suffices to consider the counterexamples 2 + 7 = 9 
and 4 + 9 = 13.) 

If P a is the largest power of a prime P that divides an integer N, i. e. P a \ N 
but P a+1 \ N, then this is denoted by P a \\N. 

We now list down several important properties of prime numbers as they 
relate to divisibility. 
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Theorem 2.1.2. If P is a prime and P \ mn, then either P \ m or P \n. 

Corollary 2.1.1. If P is a prime number and P | mim 2 • • -rrin, then P \ mi for 
some i, 1 < i < n. 

Corollary 2.1.2. // P, Qi, Q2, ■ ■ ■ , Qn are all primes and P \ Q1Q2 • • ■ Q n , then 
P = Qi for some i, 1 < i < n. 

All roads now lead to the Fundamental Theorem of Arithmetic. 
Theorem 2.1.3. Fundamental Theorem of Arithmetic 

Every positive integer N > 1 can be represented uniquely as a product of primes, 
apart from the order in which the factors occur. 

The "lexicographic representation" of a positive integer as a product of 
primes may be achieved via what is called the canonical factorization. 

Corollary 2.1.3. Any positive integer N > 1 can be written uniquely in the 
canonical factorization 

r 

N = P 1 ai P 2 a2 P r ar = Y[ P? 1 

1=1 

where, for i — 1, 2, . . . , r, each on is a positive integer and each Pi is a prime, with 
P 1 < P 2 < ... < P r . 

We illustrate these with some examples. 

Example 2.1.3. The canonical factorization of the integer 36 is 36 = 2 2 • 3 2 . 
Meanwhile, the canonical factorization for the integer 1024 is 1024 = 2 10 , while 
for 2145 it is 2145 = 3 1 • 5 1 • ll 1 • 13 1 , written simply as 2145 = 3 • 5 • 11 • 13. 
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Functions which are defined for all positive integers n are called arithmetic 
functions, or number-theoretic functions, or numerical functions. Specifically, a 
number-theoretic function f is one whose domain is the positive integers and whose 
range is a subset of the complex numbers. 

We now define three important number-theoretic functions. 

Definition 2.1.5. Let n be a positive integer. Define the number-theoretic func- 
tions d(n),a(n),(f)(n) as follows: 

d{n) = the number of positive divisors of n, 
a(n) = the sum of the positive divisors of n, 
4>{n) = the number of positive integers at most n which are also relatively prime to n. 

It would be good to illustrate with some examples. 

Example 2.1.4. Consider the positive integer n = 28. Since the positive divisors 
of 28 are 1, 2, 4, 7, 14 and 28, then by definition: 

d(28) = 6 



Note that the following list contains all the positive integers less than or equal to 
n = 28 which are also relatively prime to n: L — {1, 3, 5, 9, 11, 13, 15, 17, 19, 23, 25, 27}. 
By definition, 0(28) = 12. 



and 




d\28 
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For the first few integers, 

d(l) = 1 d(2) = 2 d(3) = 2 d(4) = 3 d(5) = 2 d(6) = 4 

while 

cr(l) = 1 <r(2) = 3 <r(3) = 4 <r(4) = 7 a(5) = 6 <r(6) = 12 

and 

0(1) = 1 0(2) = 1 0(3) = 2 0(4) = 2 0(5) = 4 0(6) = 3. 

Note that the functions d(n),cr(n), and 0(n) are not monotonic, and their 
functional values at n — 1 is also I. Also, for at least the first 3 primes p = 2, 3 
and 5, d(p) = 2, <j{p) — p + 1, and 0(p) = p — 1. 

We shall now introduce the notion of a multiplicative number-theoretic func- 
tion. 

Definition 2.1.6. A function F defined on N is said to be multiplicative if for all 

m,n EN such that gcd(m, n) = 1, we have 

F{mn) = F{m)F{n). 

Example 2.1.5. Let the function F be defined by F(n) = n k where A; is a fixed 
positive integer. Then F(mn) = {mn) k = m k n k = F(m)F(n). We have therefore 
shown that F is multiplicative. Moreover, the condition gcd(m, n) — 1 is not even 
required for the series of equalities above to hold. We call F in this example a 
totally multiplicative function. 



It turns out that the three number-theoretic functions we introduced in 
Definition 12.1.51 provide us with more examples of multiplicative functions. 



Theorem 2.1.4. The functions d, a and 4> are multiplicative functions. 



Multiplicative functions are completely determined by their values at prime 
powers. Given a positive integer n's canonical factorization 

r 

n = l[P l a , 

i=i 

then if F is a multiplicative function, we have 

r 

F(n) = ]\F(Pf<). 

This last assertion follows from the fact that prime powers derived from the canon- 
ical factorization of n are pairwise relatively prime. 



The next theorem follows from Definition 12.1.51 an d Theorem 12.1.41 as well. 

r 

Theorem 2.1.5. If n = J^Pj Qi is the canonical factorization of n > 1, then 

i=i 

r 

d(n) = JJ(oi + l), 
i=i 

i=\ i=i \ 1 / 

4>{n) = nl[(l-±). 
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We illustrate with several examples, continuing from Example 12.1.31 
Example 2.1.6. The integer 36 = 2 2 • 3 2 has 

d(36) = (2 + 1)(2 + 1) = 3 • 3 = 9 and <r(36) = (ff£)(ffr) = 7 ■ 13 = 91 
and 0(36) = 36(1 - §)(1 - §) = 12, while for the integer 1024 = 2 10 one has 

d(1024) = 10 + 1 = 11, ct(1024) = 2gzl = 2047 = 23-89 
and 0(1024) = 1024(1 - §) = 512 = 2 5 . 

Lastly, we have for the integer 2145 = 3 ■ 5 ■ 11 ■ 13 the following: 
d(2145) = (1 + 1)(1 + 1)(1 + 1)(1 + 1) = 2 4 = 16 

and 

a(2145) = (3 + 1)(5 + 1)(11 + 1)(13 + 1) = 4032 

while 

0(2145) = 2145(1 - |)(1 - |)(1 - i)(l - i) = 960. 

The following corollary follows immediately from Theorem 12.1.51 
Corollary 2.1.4. Let P be a prime number and k a fixed positive integer. Then 

d(P k ) = k + 1, a{P k ) = J2 pt = 

j=0 

and 

(f)(P k ) = P k (\ ~^) = P k -\P - 1). 

Note from Corollary 12.1.41 that for odd prime powers, the number and sum of 
divisors may or may not be prime, while <j)(P k ) is always composite for k > 1. 
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Divisibility gives rise to an equivalence relation on the set of integers, defined 
by the congruence relation. 

Definition 2.1.7. Let to be a fixed positive integer. Two integers A and B are 
said to be congruent modulo to, written as A = B (mod to), if to | (A — B); that 
is, provided that A — B = km for some integer k. When to \ (A — B), we say that 
A is incongruent to B modulo to, and we denote this by A ^ _B (mod m). 

Example 2.1.7. Let us take to = 3. We can see that 

14 = 5 (mod 3), -9 = (mod 3), and 35 = -7 (mod 3) 

because 14 - 5 = 3 • 3, -9 - = (-3) • 3, and 35 - (-7) = 14-3. 

On the other hand, 1200 ^ 2 (mod 3) because 3 does not divide 
1200 - 2 = 1198. 

We now introduce two more additional number-theoretic functions. 

Definition 2.1.8. Let n be a positive integer. Then u(n) is the number of distinct 
prime factors of n, i.e. u>(n) = Ylp^n 1 wnere eacn Pi IS prime. Furthermore, Q(n) 
is the number of primes that divide n, counting multiplicities. That is, if n has 

r 

canonical factorization n = J^-Pj a % then 

i=i 

r 

Q(n) = oti + a 2 + ■ ■ ■ + a r = a i = £*■ 

i=l Pi a ||AT 

Example 2.1.8. Let us consider n = 36 = 2 2 • 3 2 . Since it has two distinct prime 
factors (namely 2 and 3), we have ui (36) = 2. On the other hand, Q(36) = 4 since 



its total number of prime factors, counting multiplicities, is four. For m = 1024 = 
2 10 , we have w(1024) = 1 and ft(1024) = 10, while for k = 2145 = 3 • 5 • 11 • 13, 
the functions have values a; (2 145) = 4 and 11(2145) = 4. 

Notice in Example 12.1.81 that the number of distinct prime factors is less 
than or equal to the total number of prime factors (counting multiplicities). In 
general, it is true that 0(n) > u(n) for all positive integers n. 

In number theory, asymptotic density or natural density is one of the pos- 
sibilities to measure how large is a subset of the set of natural numbers N. Intu- 
itively, we feel that there are "more" odd numbers than perfect squares; however, 
the set of odd numbers is not in fact "bigger" than the set of perfect squares: 
both sets are infinite and countable and can therefore be put in one-to-one corre- 
spondence. Clearly, we need a better way to formalize our intuitive notion. 

Let A be a subset of the set of natural numbers N. If we pick randomly 
a number from the set {1, 2, ... , n}, then the probability that it belongs to A is 
the ratio of the number of elements in the set AQ{1,2, . . . ,n} and n. If this 
probability tends to some limit as n tends to infinity, then we call this limit the 
asymptotic density of A. We see that this notion can be understood as a kind of 
probability of choosing a number from the set A. Indeed, the asymptotic density 
(as well as some other types of densities) is studied in probabilistic number 
theory. 
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We formalize our definition of asymptotic density or simply density in what 
follows: 

Definition 2.1.9. A sequence a±, 02, . . . , a n with the aj positive integers and 
a,j < a J+ i for all j, has natural density or asymptotic density a, where 
< a < 1, if the proportion of natural numbers included as some aj is asymptotic 
to a. More formally, if we define the counting function A(x) as the number of a/s 
with aj < x then we require that A(x) ~ ax as x — >■ +00. 



2.2 The Abundancy Index 

As discussed in the literature review, the search for perfect numbers led 
mathematicians to consider the rather interesting quantity called the abundancy 
index/ratio. 

Definition 2.2.1. The abundancy index/ratio of a given positive integer n is 



defined as I(n) = 



cr(n 



n 



, . <r(36) 91 
Example 2.2.1. 7(36) = = — , while 

v y ' 36 36 

/(10 24) = ^i) = m 

v ' 1024 1024 

and 

a(2145) _ 4032 
2145 ~ 2145' 



1(2145) = 



We note that the abundancy index is also a multiplicative number-theoretic 
function because a is multiplicative. 
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Looking back at Definition 11.1.1), it is clear that a number N is perfect if 
and only if its abundancy index I(N) is 2. It is somewhat interesting to consider 
the cases when I(N) ^ 2. 

Definition 2.2.2. If the abundancy index I(N) < 2, then N is said to be deficient, 
while for I(N) > 2, N is said to be abundant. 

Example 2.2.2. Referring to Example 12. 2. 11 36 is abundant since 1(36) = || > 2 
while 1024 and 2145 are deficient since 7(1024) = fjjg < 2 and 7(2145) = |f§ < 2. 

Remark 2.2.1. Giardus Ruffus conjectured in 1521 that most odd numbers are 
deficient. In 1975, C. W. Anderson pQ proved that this is indeed the case by 
showing that the density of odd deficient numbers is at least ~ 0-831. On 

the other hand, Marc Deleglise (in 1998 [H]) gave the bounds 0.2474 < A{2) < 
0.2480 for the density A(2) of abundant integers. Kanold (1954) [42] showed that 
the density of odd perfect numbers is 0. 

We list down several important lemmas describing useful properties of the 
abundancy index. 



Lemma 2.2.1. 



n 




d\n 




□ 



d\n 



d\n 



d\n 



Lemma 2.2.2. If m \ n then 



a(m) 



< 



with equality occurring if and only 



m 



n 



if m = n. 
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Essentially, Lemma 12.2.21 says that any (nontrivial) multiple of a perfect 
number is abundant and every (nontrivial) divisor of a perfect number is deficient. 

Lemma 2.2.3. The abundancy index takes on arbitrarily large values. 

o~(n\) 1 n 1 

Proof. Consider the number n\. By Lemma we have — - — =y — > y —. 

n\ ^— ' a ^— ' i 

d\n\ t=l 

Since the last quantity is a partial sum of a harmonic series which diverges to 
a(n\) 

infinity, — — can be made as large as we please. □ 
n! 

Lemma 2.2.4. For any prime power P a , the following inequalities hold: 

1 P + l a(P a ) P 
1 < — ^— < „ < 



P P a P-l 

The proof of Lemma 12.2.41 follows directly from Corollary 12.1.41 

Certainly, we can find abundancy indices arbitrarily close to 1 because 
p+l 

I\P) = f° r an primes p. By Lemma T2.2.3t and since the abundancy index of 

V 

a positive integer is a rational number, one would then desire to know the "dis- 
tribution" of these ratios in the interval (1, oo). The next few results summarize 
much of what is known about the "distribution" of these ratios. 

Theorem 2.2.1. (Laatsch) The set of abundancy indices I(n) forn > 1 is dense 
in the interval (1, oo). 

However, not all of the rationals from the interval (1, oo) are abundancy 
indices of some integer. This is due to the following lemma from Weiner: 

Lemma 2.2.5. (Weiner) If gcd(m,n) = 1 and n < m < o~(n), then — is not the 
abundancy index of any integer. 
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Proof. Suppose — = — - — for some integer k. Then km = na(k) which implies 

n k 

that n | /cm, and so n \ k since m and n are coprime. Hence, by Lemma [2.2.21 we 
have 

o~(n) o-(k) m 
n ~ k n 

which yields c(n) < m - a contradiction to the initial assumption that m < 
a(n). □ 

It is now natural to define the notion of an abundancy outlaw. 

Definition 2.2.3. A rational number greater than 1 is said to be an abundancy 
outlaw if it fails to be in the range of the function I(n). 

One can use the previous lemmas to establish an equally interesting theorem 
about the distribution of abundancy outlaws. 

Theorem 2.2.2. (Weiner, Ryan) The set of abundancy outlaws is dense in the 
interval (1, oo). 

Upon inspecting the results of Theorems 12.2.11 and 12 . 2 . 2| it appears that the 
scenario for abundancy indices and outlaws is both complex and interesting. We 
shall take a closer look into the nature of abundancy outlaws in Chapter 4. 



2.3 Even Perfect Numbers 



The Greek mathematician Euclid was the first to categorize the perfect 
numbers. He noticed that the first four perfect numbers have the very specific 
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forms: 

6 = 2 1 (1 + 2) = 2- 3 
28 = 2 2 (l + 2 + 2 2 ) = 4-7 
496 = 2 4 (1 + 2 + 2 2 + 2 3 + 2 4 ) = 16 ■ 31 
8128 = 2 6 (1 + 2 + 2 2 + . . . + 2 6 ) = 64 • 127. 

Notice that the numbers 90 = 2 3 (1 + 2 + 2 2 + 2 3 ) = 8 • 15 and 
2016 = 2 5 (1 + 2 + 2 2 + . . . + 2 5 ) = 32 • 63 are missing from this list. Euclid pointed 
out that this is because 15 = 3-5 and 63 = 3 2 • 7 are both composite, whereas the 
numbers 3, 7, 31 and 127 are all prime. 

According to Book IX, proposition 36 of Euclid's Elements: Ll If as many 
numbers as we please beginning from a unit be set out continuously in double 
proportion, until the sum of all becomes a prime, and if the sum multiplied into 
the last make some number, the product will be perfect." [52J 

This observation is stated in a slightly more compact form as follows: 

Theorem 2.3.1. (Euclid) If 2 n — 1 is prime, then N = 2 n_1 (2 n — 1) is perfect. 

Proof. Clearly the only prime factors of N are 2 n — 1 and 2. Since 2 n — 1 occurs 
as a single prime, we have simply that o~(2 n — 1) = 1 + (2 n — 1) = 2 n , and thus 

a(N) = a(2 n - 1 )a(2 n - 1) = (^)2 n = 2 n (2 n - 1) = 2N . 

Therefore, N is perfect. □ 



The task of finding perfect numbers, then, is intimately linked with finding 
primes of the form 2 n — 1. Such numbers are referred to as Mersenne primes, after 
the 17th-century monk Marin Mersenne, a contemporary of Descartes, Fermat, 
and Pascal. He investigated these unique primes as early as 1644. Mersenne knew 
that 2" — 1 is prime for n — 2, 3, 5, 11, 13, 17, and 19 - and, more brilliantly, 
conjectured the cases n = 31, 67, 127, 257. It took almost two hundred years to 
test these numbers. 

There is one important criterion used to determine the primality of Mersenne 
numbers: 

Lemma 2.3.1. (Cataldi- Fermat) If 2 n — 1 is prime, then n itself is prime. 

Proof. Consider the factorization of x n — 1 = (x — l)(x n ~ l + . . . + x + 1). Suppose 
n = rs, where r, s > 1. Then 2 n - 1 = (2 r ) s - 1 = (2 r - l)((2 r ) s - 1 + . . . + 2 r + 1), 
so that (2 r — 1) | (2 n — 1) which is prime, a contradiction. □ 

Note that the converse of Lemma f2.3. II is not true - the number 2 11 — 1 which 
is equal to 2047 = 23-89 is composite, yet 11 is prime, for instance. 

Should all perfect numbers be of Euclid's type? Leonard Euler, in a posthu- 
mous paper, proved that every even perfect number is of this type. [48] 

Theorem 2.3.2. (Euler) If N is an even perfect number, then N can be written 
in the form N = 2 n_1 (2™ — 1), where 2 n — 1 is prime. 

Proof. Let iV = 2 n ~ 1 m be perfect, where m is odd; since 2 does not divide m, it 
is relatively prime to 2 n_1 , and 



a(N) = a{2 n - l m) = a{2"- l )a{m) = (^)cx(m) = (2 n - l)a{m). 

N is perfect so cr(N) = 2N = 2(2"~ 1 m) = 2™m, and with the above, 
2 n m = (2™— l)a(m). Since 2 n — 1 is odd, (2 n — 1) | m, so we can write m = (2 n —l)k. 
Now (2"-l)(x(m) = 2"(2 n -l)A;, which implies a{m) = 2 n k = (2 n -l)k+k = m+k. 
But k | m so <j(m) = m + k means m has only two (2) divisors, which further 
implies that k = 1. Therefore, a(m) = m + 1 and m is prime. Since (2™ — l)\m, 
2 n — 1 = m. Consequently, iV = 2 n_1 (2 n — 1) where 2 n — 1 is prime. □ 

Even perfect numbers have a number of nice little properties. We list down 
several of them here, and state them without proof [71] : 

• If iV is an even perfect number, then N is triangular. 

• UN — 2 n ~\2 n - 1) is perfect then iV = l 3 + 3 3 + . . . + (2^ - l) 3 . 

• If N = 2 n ~ 1 (2 n — 1) is perfect and N is written in base 2, then it has 2n — 1 
digits, the first n of which are unity and the last n — 1 are zero. 

• Every even perfect number ends in either 6 or 8. 

• ( Wantzel) The iterative sum of the digits (i.e. digital root) of an even perfect 
number (other than 6) is one. 

Today 44 perfect numbers are known, 2 88 (2 89 — 1) being the last to be 
discovered by hand calculations in 1911 (although not the largest found by hand 
calculations), all others being found using a computer. In fact computers have led 
to a revival of interest in the discovery of Mersenne primes, and therefore of perfect 
numbers. At the moment the largest known Mersenne prime is 2 32582657 — 1. It was 



discovered in September of 2006 and this, the 44th such prime to be discovered, 
contains more than 9.8 million digits. Worth noting is the fact that although this 
is the 44th to be discovered, it may not correspond to the 44th perfect number as 
not all smaller cases have been ruled out. 

2.4 Odd Perfect Numbers 

The Euclid-Euler theorem from Section 2.3 takes care of the even perfect 
numbers. What about the odd perfect numbers? 

Euler also tried to make some headway on the problem of whether odd 
perfect numbers existed. He proved that any odd perfect number N had to have 
the form 

N={4m + l) 4k+1 b 2 
where 4m + 1 is prime and gcd(4m + 1,6) = 1. 

In Section 2.3, we have followed some of the progress of finding even perfect 
numbers but there were also attempts to show that an odd perfect number could 
not exist. The main thrust of progress here has been to show the minimum number 
of distinct prime factors that an odd perfect number must have. Sylvester worked 
on this problem and wrote: 



. . . the existence of [an OPN] - its escape, so to say, from the complex 
web of conditions which hem it in on all sides - would be little short of a miracle. 



(The reader is referred to Section 1.1.2 of this thesis for a survey of the most 
recent conditions which an OPN must satisfy, if any exists.) 

We give a proof of Euler's characterization of OPNs here |48j : 

Theorem 2.4.1. (Euler) Let N be an OPN. Then the prime factorization of N 
takes the form N = q +1 P i 2ai . . -p r 2ar ; where q = 1 (mod 4). 

Proof. Let N = l\ ei l 2 e2 • ■ ■ l s es for some primes li, l 2 , ■ ■ ■ , l s - Since N is odd, all 
are odd. Finally, a(N) = 2N. Since a(N) = a(h e H 2 e2 ■ ■ ■ l s €s ) = o-{h ei )a{l 2 e2 ) ■ ■ ■ a(l 
we take a look at cr(/ e ) = 1 + I + I 2 + . . . + l e , a sum of e + 1 odd numbers. 
This is odd only if e is even. Since a(h ei l 2 e2 ■ ■ ■ l s e °) = cx(h ei )a(l 2 e2 ) ■ ■ ■ <x(//*) = 
2l 1 ei l 2 e2 ■ ■ ■ l s e % we can only get one factor of 2. So the are even, all except one, 
say e x . So iV = l 1 ei Pl 2a ^ ■ --p 2 ^. 

We have 2 | cx(Zi ei ) but 4 \ a{li ei ). Since l\ is odd, e\ is odd. Now, modulo 
4, we see that either l x = 1 (mod 4) or l x = — 1 (mod 4). But if Zx = — 1 (mod 4), 
then 

a(h ei ) = l + h + h 2 + h 3 + ... + l^- 1 + h ei 
= 1 + (-1) + 1 + (-1) + . . . + 1 + (-1) = (mod 4), 

which is clearly a contradiction since 4 { <j(Zi ei ). Thus, l\ = 1 (mod 4). Now, 
a(h ei ) = 1 + h + h 2 + h 3 + ... + h ei - 1 + h ei = 1 + 1 + 1 + 1 + .. . + 1 + 
1 = ei + 1 (mod 4). Since e\ is odd, either e\ + 1 = (mod 4) or e\ + 1 = 2 
(mod 4). If ei + 1 = (mod 4), then 4 | a{l\ ei ) which is again a contradiction. 



So t\ + 1 = 2 (mod 4) <^ ei + 1 = 4e + 2, that is, t\ = 4e + 1. Consequently, 

jy = q ie+l p 2 ai . . .p r 2a r; f or g = 1 ( mQc l 4). □ 

We call q in Theorem 12.4. II the special/Euler prime of N, while g 4e+1 will 
be called the Euler's factor of N. 

Interestingly, it is possible to show that no two consecutive integers can be 
both perfect, using the Euclid- Euler theorem on the form of even perfect numbers 
from Section 2.3 and Euler's characterization of odd perfect numbers in this sec- 
tion. We shall give a discussion of the proof of this interesting result in Chapter 
3. Moreover, we shall give there the (easy) proof of the fact that an odd perfect 
number must be a sum of two squares. We also give congruence conditions for 
the existence of odd perfect numbers in the next chapter. For the most part, the 
proofs will be elementary, requiring only an intermediate grasp of college algebra 
and the concepts introduced in this chapter. 



Chapter 3 

OPN Solution Attempts 1: 
Some Old Approaches 

In the previous chapters, we derived the possible forms of even and odd perfect 
numbers, and also surveyed the most recent results on the conditions necessitated 
by the existence of OPNs. In this chapter, we introduce the reader to the fla- 
vor of the mathematical techniques used to formulate theorems about OPNs by 
researchers who lived prior to the 21st century. 

The following are some of the traditional attempts made by mathematicians 
(both amateur and professional) to prove or disprove the OPN Conjecture in the 
pre-21st century: 

• Increasing the lower bound for the number of distinct prime factors, ou(N), 
that an OPN N must have; 

• Increasing the lower bound for the magnitude of the smallest possible OPN, 
if one exists; 

• Deriving congruence conditions for the existence of OPNs. 
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All these itemized approaches attempt to derive a contradiction amongst 
the stringent conditions that an OPN must satisfy. 



3.1 Increasing the Lower Bound for uj(N) 



Recall that, from Lemma I2.2.4[ we have the following strict inequality for the 
abundancy index of a prime power: — - < ' 



pa p — 1 
This gives rise to the following lemma: 

oj(N) 

Lemma 3.1.1. If N is a perfect number with canonical factorization N = J^J P% ai , 

i=i 

then 

lu(N) n u(N) 

2 < 



II/: Hi rSi 



i=l " i=l 

Proof. This follows as an immediate consequence of Lemma 12.2.41 and the defini- 
tion of perfect numbers. □ 

Note that Lemma 13.1.11 applies to both even and odd perfect numbers. 

The following is another very useful lemma: 

Lemma 3.1.2. If N is a perfect number with canonical factorization N = |~J , 

i=i 

then 



i=l 

where < fa < a^Vi. 



1 



Pi " ' B ft 
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Proof. This immediately follows from the definition of perfect numbers and the 
fact that the abundancy index for prime powers is an increasing function of the 
exponents. □ 

When considering OPNs, Lemmas 13.1.11 and 13.1.21 are very useful because 
they can yield lower bounds for u(N). Indeed, it was in using these lemmas 
(together with some ingenuity) that pre-21st century mathematicians were able 
to successfully obtain ever-increasing lower bounds for the number of distinct 
prime factors of an OPN N. 

We now prove the classical result: u An OPN must have at least three distinct 
prime factors" . 

Theorem 3.1.1. If N is an OPN, then u(N)>3. 

Proof. Let N be an OPN. Since prime powers are deficient, u(N) > 2. Suppose 
u(N) = 2. Since N is odd, P x > 3 and P 2 > 5 where N = P 1 ai P 2 a2 is the 
canonical factorization of N. Using Lemma 13.1.11 

"(AT) 2 

i=i i=i 
Pi \ ( Pi \ i i ii 



p x -\) \p 2 -ij i-^i-^-i-ii-i 

11 3 5 15 
= II=24 = T = 1 - 875 <2 

3 5 

Thus, the assumption oj(N) = 2 for an OPN iV has resulted to the contra- 
diction 2 < 2. This contradiction shows that u(N) > 3. □ 
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More work is required to improve the result of Theorem 13. 1 . II to u(N) > 4, 
if we are to use a similar method. 

Theorem 3.1.2. If N is an OPN, then u(N) > 4. 

Proof. Let N be an OPN. By Theorem EHH u(N) > 3. Assume u(N) = 3. Since 
N is odd, Pi > 3, P 2 > 5, and P 3 > 7 where N = p^p^p^ i s the canonical 
factorization of N. Using Lemma [3. 1.11 

"(AT) 3 

^ II,. . II/. 

8=1 1 1=1 1 

Pt \ ( P 2 \ ( P 3 \ 1 1 1 < 1 1 1 



Pi-y\P2-iJ\Pz-v i-^i-^i-^ - i-ii-ji-i 

= Ui = -~ = - = - = 2.1875 
HI 246 48 16 

3 5 7 

whence we do not arrive at a contradiction. Now, suppose P\ > 5, P 2 > 7, and 
P 3 > 11. Using Lemma 13.1.11 again, we get 2 < (|)(|)(^) = ||, which is a 
contradiction. Consequently, we know that Pi = 3. Then, with this additional 
information about Pi, if we assume that P 2 > 7, we arrive at 2 < (|)(|)(y§) = 
f§5 = 15 use °f Lemma 13. 1 . 1[) which is again a contradiction. Hence, we also 
know that P 2 = 5. Thus, P3 > 7. Furthermore, using Lemma 13.1.11 again, the 
following inequality must be true: 



2 4 P, - 1 15 P 



3 



Solving this last inequality, we get P 3 < 16. This inequality, together with the 
fact that P 3 is prime, gives us 3 possible cases to consider: 
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Case 1: N = 3 ai 5 a2 7° 3 

Since N is odd, 4 cannot divide a(N) = a(3 ai )a(5 a2 )a(7 Q3 ) = 2N = 2-3 ai 5 a ' 2 7 as . 
In particular, 4 \ a(3 ai ) and 4 \ a(7 a[i ). But <x(3 Ql ) = 4 for e*i = 1 and a(7 as ) = 8 
for a 3 = 1. Consequently, «i > 2 and a 3 > 2. Now, by using Lemma [3.1.21 with 
f3\ = 2, (3 2 = 1 and /3 3 = 2, we have: 



1 1 \ / 1\ / 1 1 \ 494 

2> 1 + - + 1 + T 1 + - + 



3 3VV 57V 7 7V 245' 

which is a contradiction. Hence, there is no OPN of the form iV = 3 Q i5 a 2 7 a 3. 
Case 2: N = 3 Q i5 a2 ll Q3 

„ , , , 3(r(5 a2 )ll 3 611 99 , , 

^4. Using Lemma 3.1.1 with «o = 1 gives 2 < = = — , which is 

6 6 2 5^ 10 2510 50 

a contradiction. 

B. Let a 2 > 2. For at = 1, 2, 3, we have a{3 ai ) = 4, 13, 40. Since 4 f <r(JV), a x ^ 1 
and ai ^ 3. Also, 13 f a(N) = 2N = 2 • 3 ai 5 a2 ll a3 and thus, 13 \ cr(3 ai ). Hence, 
ai 7^ 2, which implies that oti > 4. Similarly, <r(ll a3 ) = 12 for a 3 = 1, which 
contradicts the fact that 4 \ cr(ll a3 ). Thus, a 3 > 2. Now, using Lemma [3.1.21 with 

fa = 4, 02 = 2 and /3 3 = 2, we get: 



1 1 1 1 \ / 1 1 \ / 1 1\ 4123 
2> l + - + ^ + ^ + ^r l + - + -7 l + — + -o 



3 3 2 3 3 3 4 J \ 5 5 2 J V 11 H 2 / 2025 



which is a contradiction. Hence, there is no OPN of the form iV = 3 ai 5 a2 ll a3 . 



Case 3: N = 3 ai 5 a2 13 a3 

. TT . T . , , . n 3(t(5 Q2 ) 13 3 6 13 39 . 

A. Using Lemma 3.1.1 with a 2 = 1 gives 2 < = = — , which is 

6 6 2 5« 2 12 2 5 12 20 

a contradiction. 
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B. Let a 2 > 2. Similar to what we got from Case 2B, we have a\ ^ 1, 3. 

1. ax = 2 

, , , ct(3 Qi ) 5 13 13 5 13 845 , , 

Using Lemma 3.1.11 we have: 2 < = = , which is a con- 

6 1 3 Q i 412 9 412 432' 

tradiction. 

2. q;i > 4 

Suppose a 3 = 1. Then cx(13 a3 ) = 14 = 2 • 7, yet 7 f tr(JV) = 2iV = 2 ■ 3 Ql 5 Q2 13 a3 . 
Therefore, 7{ a^" 3 ), and then 03 > 2. Now, by using Lemma 13. 1.21 with (3\ = 4, 

/3 2 = 2 and /3 3 = 2, we get: 



1 1 1 1 \ / 1 1 \ / 1 I \ 228811 
2> l + - + ^ + rr + -r l + - + -r l + — + -0 



3 3 2 3 3 3 4 y V 5 5 2 y V 13 132 / H4075' 
which is again a contradiction. Hence, there is no OPN of the form N = 3 ai 5" 2 13" 3 . 

Therefore, there is no OPN with exactly three (3) distinct prime factors, i.e. 
an OPN must have at least four (4) distinct prime factors. □ 

Remark 3.1.1. Using more recent findings on an upper bound for OPNs by 

Nielsen (50] and on a lower bound by Brent et al. [4], it is possible to extend the 

results in this section. Thus, 10 300 < N < 2 4 " (iV) , and this gives 

> p + log(3) - log(log(2m > 
V log(4) J 

which implies that u(N) > 5 since u(N) should be an integer. (There is a project 



underway at http:/ /www.oddperfect.org, organized by William Lipp, which hopes 
to extend the lower bound for OPNs to 10 500 , or 10 600 even.) Indeed, even more re- 
cently (2006), Nielsen [51] was able to show that u(N) > 9, "ultimately [avoiding] 
previous computational results for [OPNs]". 
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3.2 Increasing the Lower Bound for an OPN N 

We begin with a very useful inequality (which we shall not prove here) that can 
yield our desired estimates for a lower bound on OPNs: 

Lemma 3.2.1. The Arithmetic Mean-Geometric Mean Inequality 

Let {Xi} be a sequence of (not necessarily distinct) positive real numbers. Then 

the following inequality must be true: 



n \ n 



i=l \i=l 



Equality holds if and only if all of the Xi 's are equal. 



We now derive a crude lower bound for an OPN N = j^J P^ in terms of 



the cKj's: 



Lemma 3.2.2. Let A ] [ Pf* be an OPN. Then 



i=i 



n («« + i ) 



N > 



i=l 



V 



/ 



Proof. a{P i ° li ) = VjPj- 7 . Applying Lemma [3.2.11 and noting that prime powers 



3=0 

of the P^s are distinct, we have for each i: 

j 



^— ' a- + 1 

5]jy>(a i + l)nfl a * + 1 =(ai + lW ==0 1 



i=o 



i=o 
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= («, + 1)P, + X ) = (a 4 + l)pT. 
Consequently, by multiplying across all i: 

w{N) ui(N) Ot-i ui{N) OCi 

I [ a{Pi ai ) > JJ (a 4 + 2 . But V^V = JJ Pi 2 , which means that 

i=l i=l i=l 

n<w) n (a* + 1) n ^ n^+i) 

N = — > — — = VN- l=1 



2 2 2 

Solving this last inequality for iV gives us the desired result. □ 

w(JV) 

Remark 3.2.1. In the canonical factorization = J^J p Qi of an OPN N, since 

i=i 

ot-i > 1 for all z, from Lemma [3.2.21 we have the crude lower bound iV > 2 2aj ( 7V )~ 2 , 
which, together with Nielsen's lower bound of u(N) > 9 for the number of distinct 
prime factors of N, yields the lower bound iV > 2 16 = 65536 for the magnitude of 
the smallest possible OPN. This lower bound can, of course, be improved. Indeed, 
Brent, et al. [1] in 1991 showed that it must be the case that iV > 10 300 . 

Remark 3.2.2. Note that nowhere in the proof of Lemma [3.2.21 did we use the 
fact that iV is odd. Hence, Lemma 13.2.21 applies to even perfect numbers as well. 

We can make use of the results on the lower bound for the number of distinct 
prime factors (latest result is at 9 by Nielsen) , lower bound for the smallest prime 
factor (currently at 3 - mathematicians are still unable to rule out the possibility 
that an OPN may be divisible by 3), and the nature of the exponents (the spe- 
cial/Euler prime has the sole odd exponent while the rest of the primes have even 
exponents) to derive a larger lower bound for an OPN. 



43 



ui(N) 



Since all, except for one, of the distinct prime factors of an OPN N 



I [ Pf* have even exponents, then a, > 2 for all but one z, say i = j, for which 

i=l 

ctj > 1. (Note that <x,- is the exponent of the special or Euler prime.) Thereupon, 
we have the following improvements to the results in Remark 13.2.11 

n + !) > 3 " (jvki ■ 2 
n («* + r 



A/" > 



V 



> 3 



2w(JV)-2 



/ 



w(iV) > 9 [Nielsen] N > 3 16 = 43046721 



Note that there is approximately a 655.84% improvement in the magnitude of the 
bound thus obtained for the smallest possible OPN as compared to the previous 
crude lower bound of 2 16 . The novelty of the approach of Lemma 13.2.21 can be 
realized if we consider the fact that we did NOT need to check any of the odd 
numbers below 43046721 to see if they could be perfect. 

Remark 3.2.3. We casually remark that the lower bound of 3 16 obtained for an 
OPN here improves on the classical bound of 2 • 10 6 obtained by Turcaninov in 
1908. However, at that time, the best-known bound for the number of distinct 
prime factors of an OPN iV was u(N) > 5, which was shown to be true by 
Sylvester in 1888, whereas we used the bound u(N) > 9 by Nielsen (2006) here. 



As mentioned in Remark 13.2.31 A. Turcaninov showed in 1908 that no odd 
number less than two million can be perfect. The figure 2- 10 6 is generally accepted 
as the minimum in standard texts. Nonetheless, it is easy to show by means of 



well-known proofs that the smallest possible OPN must be greater than ten billion 
(i.e. 10 10 ). 

Let the prime factorization of an OPN N be given by 

N = Pi ai Qi bl Q2 b2 ---Q m bm 

where a\ is odd and 61, 62, • • • , b m are even. The following conditions must hold: 

• Euler had shown that P x = ai = 1 (mod 4). 

• Sylvester demonstrated that it must be the case that m > 4. 

• Steuerwald proved that b\ = b 2 = ■ ■ ■ = b m = 2 is not possible. 

• Brauer extended the last result to b; L 7^ 4, when b\ — b 2 — • • • — fej-i = 
bi+i — ■ ■ ■ — b m — 2. 

• Sylvester also showed that 105 = 3-5-7 does not divide A^ . (We give a 
proof of this result in Section 3.4.) 

The only numbers less than ten billion which satisfy all these itemized con- 
ditions are 3 6 • 5 2 • ll 2 • 13 2 • 17 and 3 6 • 5 2 • ll 2 • 13 • 17 2 . One can verify that each 
of these two is abundant by directly computing the sum of its divisors. 

We end this section with a copy of an email correspondence between the 
author and Richard Brent, one of the three co-authors of the 1991 paper which 
showed that an OPN N must be bigger than 10 300 . The significance of the email 
lies with the fact that the author of this thesis was able to show that N > p k a(p k ) 



where k is unrestricted (i.e. k can be even or odd). In fact, the author was able 
to show the slightly stronger statement N > ^p k a(p k ) with k unrestricted. (Al- 
though Dr. Brent did not mention it in his reply, this result could give a higher 
lower bound for OPNs.) We shall give a proof of this last result in Theorem 4.2.5. 
For now, let us take a look at the author's email inquiry and how Dr. Brent 
responded to it: 
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From: 



Richard Brent [rpbrent@gmail.com] 



Sent: 



Tuesday, December 04, 2007 9:38 AM 



To: 



Jose Dris 



Subject: Re: Inquiry regarding the lower bound of 10 that you obtained for 
odd perfect numbers 

Dear Arnie, 

On 22/11/2007, Jose Dris (Jose.Dris@safeway.com) wrote: 
Hi Dr. Brent, 

I am Arnie Dris, a candidate for the degree of MS in Mathematics at DLSU 
here at Manila, Philippines. I am currently in the process of writing up a thesis on 
odd perfect numbers, and I came across your lower bound of 10 300 for odd perfect 
numbers that you obtained together with two co-authors. 

In your paper, you stated that you used the simple observation that 
N > p k a(p k ) where p k \\N and k is even. I would just like to ask whether this 
observation was used in ALL cases that you have considered, thereby proving that 
N > 10 300 in each case? 



It was a long time ago, but as I recall we used that observation in most cases. 
There were a few "hard cases" where we could not compute the sigma function 
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a(p k ) because we did not know the complete factorisation that is needed to do 
this, e.g. cr(3221 42 ) was a 148-digit composite number that we could not factorise 
at the time (it may have been factored since then). In such cases we had to use 
a more complicated method. See the "Proof of Theorem 1" in the paper. 

You can get the paper online at 

http:/ /wwwmaths. anu.edu.au/~brent/pub/publl6. html 

and there's also a link there to the computer-generated "proof tree". 

William Lipp has a project to extend the bound 10 300 by much the same method 
but with more factorisations (since computers are faster now and algorithms such 
as the number field sieve are available). He hopes to go at least to 10 400 and 
possibly further. 

Regards, 
Richard Brent 

Prof R. P. Brent, ARC Federation Fellow 
MSI, ANU, Canberra, ACT 0200, Australia 
rpbrent@gmail.com 

|http: / /wwwmaths .anu.edu.anu/ ~brent/ 
phone: +61-4-18104021 
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3.3 Congruence Conditions for an OPN N 



In this section, we revisit a theorem of Jacques Touchard dating back from 1953. 
Touchard proved that any odd perfect number must have the form 12m + 1 or 
36m + 9. His proof relied on the fact that the numbers a(k) satisfy 



a recursion relation derived by Balth. van der Pol in 1951 using a nonlinear partial 
differential equation. We give here Judy Holdener's proof of the same result in 
2002, which is much shorter and more elementary than Touchard's proof. The 
proof was inspired by Francis Coghlan's solution to Problem #10711 published in 
the American Mathematical Monthly in 2001 regarding the nonexistence of two 
consecutive perfect numbers. 

First, we briefly spell out a lemma on a congruence condition for an OPN 



Lemma 3.3.1. If N = 5 (mod 6), then N is not perfect. 

Proof. Assume that N = 5 (mod 6). Then N is of the form Qk + 5 = 3(2k + 1) +2, 
so N = 2 (mod 3). Since all squares are congruent to 1 modulo 3, iV is not a 




N: 



square. Further, note that for any divisor d of N, N = d ■ (^) = 2 = —1 (mod 3) 
implies that either d = — 1 (mod 3) and ^ = 1 (mod 3), or d = 1 (mod 3) and 
f = -1 (mod 3). Either way, d + (f ) = (mod 3), and 




d\N, d<VN 



d+ — \ =0 (mod 3). 
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Therefore, cr(N) = (mod 3) while 2N = 4=1 (mod 3). These computations 



Using a similar argument, we can also show Euler's result that any OPN 
is congruent to 1 modulo 4. For suppose otherwise that N = 3 (mod 4). Then, 
again, N is not a square and 



cr(N) = frf+^-j =0 (mod4). 

d\N, d<^N 

Hence, a(N) = (mod 4), while 2N = 6 = 2 (mod 4). 

Lemma 13.3.11 generalizes immediately to the following Corollary: 

Corollary 3.3.1. If M is a number satisfying M = 2 (mod 3), i/ien M zs noi 
perfect. 

Proof. Note that Lemma 13.3.11 takes care of the case when M is odd. We now 
show that the statement is true for even M. 

To this end, suppose M = 2 (mod 3) is even. We show here that M cannot 
be perfect. Suppose to the contrary that M is even perfect. Then by the Euclid- 
Euler Theorem, M = 2 P_1 (2 P — 1) for some prime number p. If p = 2, then 
M = 6 which is divisible by 3. Assume p > 3. Then p is an odd prime. Thus, 
2 p = (-l) p = (-l) p -\-l) = 1 • (-1) = -1 = 2 (mod 3). Therefore, 2 P - 1 = 1 
(mod 3). Also, 2 P_1 = ( — l) p_1 = 1 (mod 3) since p is assumed to be an odd 



show that N = Qk + 5 cannot be perfect. 



□ 




prime. 



Thus, if M is to be even perfect, either M = (mod 3) (which occurs 
only when p = 2 and M = 6) or M = 1 (mod 3) (when p is a prime > 3). 
Consequently, if M = 2 (mod 3) is even, it cannot be perfect. Lemma [3.3.11 says 
that this is also true when M is odd. Hence, we have the general result: If a 
number M satisfies M = 2 (mod 3), then M cannot be perfect. □ 

We we will use the following formulation of the Chinese Remainder Theorem 
to prove the next major result: 

Theorem 3.3.1. Chinese Remainder Theorem 

Suppose ni,n 2 , . . . ,Uf, are integers which are pairwise relatively prime (or co- 
prime). Then, for any given integers ai, a%, . . . , a^, there exists an integer x solving 
the system of simultaneous congruences 

x = a\ (mod n\) 
x = a 2 (mod n 2 ) 

x = a k (mod n^) 

Furthermore, all solutions x to this system are congruent modulo the product n = 
n\ri2- • -Tik- Hence x = y (mod rii) for all 1 < % < k, if and only if x = y 
(mod n). 

Sometimes, the simultaneous congruences can be solved even if the rii 's are 
not pairwise coprime. A solution x exists if and only if = aj (mod gcd(nj, rij)) 
for all % and j. All solutions x are then congruent modulo the least common 
multiple of the n^. 



We can now use Lemma 13.3.11 to prove Touchard's theorem. 

Theorem 3.3.2. (Touchard) An OPN must have the form 12m + 1 or 36m + 9. 

Proof. Let iV be an OPN and apply Lemma [3.3.11 Any number of the form 6k + 5 
cannot be perfect, so N must be of the form 6k + 1 or 6k + 3. But from Euler's 
result, we know that N is of the form Aj + 1. Hence either N = 6k + 1 and 
iV = Aj + 1, or iV = 6k + 3 and N = Aj + 1. We now attempt to solve these two 
sets of simultaneous equations for N, thereby deriving congruence conditions for 
N: 

Case 1: N = 6k + 1 and N = A j + 1. This means that iV - 1 = 6k = Aj = 
LCM(4, 6)m = 12m (by the Chinese Remainder Theorem) where k = 2m and 
j = 3m, which implies that iV = 12m + 1. 

Case 2: iV = 6k + 3 and N = Aj + 1. This means that N + 3 = 6k + 6 = 
Aj + 4 = 6(k + 1) = A(j + 1) = LCM(A, 6)m = Up (by the Chinese Remainder 
Theorem) where k + 1 = 2p and j + 1 — 3p, which implies that N = 12p — 3. On 
setting p = m + 1, we get iV = 12m + 9. 

Finally, note that in Case 2, if N — 12mo + 9 and 3 \ mo, then a(N) = 
(j(3(4m + 3)) = cr(3)cr(4m + 3) = 4or(4m + 3). With this, we have a(N) = 
(mod 4), while 2N = 2(12m + 9) = 24m + 18 = 4(6m + 4) + 2 = 2 (mod 4). 
Therefore, iV cannot be perfect if 3 f m in Case 2, and we conclude that 3 | mo in 
this case, and on setting m = 3m, we get iV = 12m + 9 = 12(3m) + 9 = 36m + 9. 

□ 



Remark 3.3.1. We emphasize that Touchard's theorem is really simple. Hold- 
ener's proof as presented here is indeed elementary; it does not make use of the 
concept of unique factorization nor of sigma multiplicativity (other than in show- 
ing that 3 divides m when 12m + 9 is perfect). Touchard's result emerges after 
summing divisors in pairs, and this can always be done because perfect numbers 
are never squares. 

In January of 2008, Tim Roberts made a post at 

http:// www, unsolvedproblems . org/ UP/ Solutions . htm | 
where he outlined an improvement to Theorem 13. 3. 2[ 

Theorem 3.3.3. (Roberts) Let N be an OPN. Then either one of the following 
three congruences must hold: 

• N = 1 (mod 12). 

• N = 117 (mod 468). 

• TV = 81 (mod 324). 

Proof. Let N be an OPN. We note that, if 3 | N, then 3 fe | N, where k is even 
(Euler). If k — 0, then by Theorem I3.3.2[ N = 1 (mod 12). Also, by the factor 
chain approach, if N is an OPN and a factor of N is 3 k , then N is also divisible 
by a(3 fe ) = 1 + 3 + 3 2 + . . . + 3 fc . If k = 2, then again by Theorem I31T21 N = 9 
(mod 36). Further, since N is an OPN, we know that cr(3 2 ) = 1 + 3 + 3 2 = 13 
divides N. Hence, N = (mod 13). From the Chinese Remainder Theorem, we 
can deduce that N = 117 (mod 468). If k > 2, then N is divisible by 3 4 = 81. 



Thus, (again by Theorem I3.3.2p iV must satisfy both N = 9 (mod 36) and N = 
(mod 81). Again, from the Chinese Remainder Theorem, we can deduce that 
N = 81 (mod 324). Thus, if N is an OPN, then either N = 1 (mod 12), JV = 117 
(mod 468) or iV = 81 (mod 324). □ 

It is, of course, similarly possible to further refine the last of these results, 
by separately considering even values of k bigger than 4. 

3.4 Some Interesting Results on Perfect Num- 
bers 

We conclude this chapter with the following (interesting) results on perfect num- 
bers (with emphasis on OPNs): 

• No two consecutive integers can be both perfect. 

• An odd perfect number cannot be divisible by 105. 

• An odd perfect number must be a sum of two squares. 

3.4.1 Nonexistence of Consecutive Perfect Numbers 

From Corollary 13.3.11 we see that a number M (odd or even) satisfying 
M = 2 (mod 3) cannot be perfect. We make use of this observation to prove that 
no two consecutive integers can be both perfect. (This was shown to be true by 
James Riggs and Judy Holdener through a joint undergraduate research project 
in 1998.) 



First, suppose N is an OPN. By Euler's characterization of OPNs, N = 1 
(mod 4). We claim that N + 1 cannot be an even perfect number. Observe that 
N + l-2 (mod 4) means that 2 | (N + 1) but 4 f (JV + 1). The only even perfect 
number of the form N + 1 = 2 P ~ 1 (2 P — 1) satisfying these two conditions is the 
one for p = 2, i.e. N + 1 = 6. But this means that, by assumption, N = 5 must 
be an OPN, contradicting the fact that N = 5 is deficient. 

Next, we also claim that N — 1 cannot be (even) perfect, if N is an OPN, 
where N — 1 = (mod 4). From the discussion of the proof of Corollary 13.3.11 
since N — 1 = (mod 4) it follows that N — 1 = 2 P ~ 1 (2 P — 1) for some primes 
p and 2 P — 1 with p > 3. Thus, N — 1 = 1 (mod 3), which implies that N = 2 
(mod 3). But our original assumption was that iV is an OPN, contradicting the 
criterion in Corollary 13.3.11 Consequently, this means that iV — 1 is not perfect 
in this case. 

We have shown in the preceding paragraphs that, if N is an OPN of the 
form 4m + 1, then it cannot be true that N — 1 or iV + 1 are also (even) perfect. 
To fully prove the assertion in the title of this section, we need to show that iV — 1 
and N + 1 cannot be OPNs if iV is an even perfect number. 

To this end, suppose N is even perfect, that is, iV = 2 P_1 (2 P — 1) for some 
primes p and 2 P — 1. If p = 2, then N = 6, and clearly, N — 1 = 5 and N + 1 = 7 
are not perfect since they are both primes (and are therefore deficient). 



Now let p be a prime which is at least 3. Then N = (mod 4), whence 
it follows that N — 1 = 3 (mod 4) and N — 1 cannot be an OPN by Euler's 
characterization. Also, from the proof of Corollary 13.3. 1[ note that if N is an even 
perfect number with p > 3, then N = 1 (mod 3), which implies that iV + 1 = 2 
(mod 3), which further means that N + 1 cannot be perfect by the criterion in 
Corollary [3XU 

3.4.2 OPNs are Not Divisible by 105 

Mathematicians have been unable, so far, to eliminate the possibility that 
an odd perfect number may be divisible by 3. However, by use of Lemma I3.1.2[ 
we can show that an odd number divisible by 3, 5 and 7 cannot be perfect. (This 
was proved by Sylvester in 1888.) 

To this end, suppose that N is an OPN that is divisible by 3, 5, and 7. 
Then iV takes the form N = 3 a 5 b 7 c n. Suppose that 5 is the special/Euler prime 
of N, so that b > 1. By Euler's characterization of OPNs, a and c must be even, 
so we may take a > 2 and c > 2. Without loss of generality, since u(N) > 9 
(see Remark I3.1.ip we may safely assume that n > ll 2 . Then J(n) > 1, and 
consequently, by use of Lemma 13.1.2^ we have: 

13 6 57 4446 

> 1 = > 2.0163 > 2 

9 5 49 2205 

resulting in the contradiction 2 > 2. 



If, in turn, we assume that 5 is not the special/Euler prime of N (so that 
b > 2), then without loss of generality, since n must contain the special/Euler 
prime and u(N) > 9, we can assume that n > 13. In this case, it is still true that 
I{n) > 1 (also that a > 2 and c > 2). Hence, by use of Lemma 13.1.21 we obtain: 

resulting, again, in the contradiction 2 > 2. 

We are therefore led to conclude that an OPN cannot be divisible by 3-5-7 = 

105. 

3.4.3 OPNs as Sums of Two Squares 

We borrow heavily the following preliminary material from 
http: / / en.wikipedia.org/ wiki/Proofs_of_Fermat's_theorem_on_sums_of_two_squares 
This is in view of the fact that the special/Euler prime p of an OPN N = p k m 2 
satisfies p = 1 (mod 4). 

Fermat's theorem on sums of two squares states that an odd prime p can 
be expressed as p = x 2 + y 2 with x and y integers if and only if p = 1 (mod 4). 
It was originally announced by Fermat in 1640, but he gave no proof. The only 
if clause is trivial: the squares modulo 4 are and 1, so x 2 + y 2 is congruent to 
0, 1, or 2 modulo 4. Since p is assumed to be odd, this means that it must be 
congruent to 1 modulo 4. 
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Euler succeeded in proving Fermat's theorem on sums of two squares in 
1747, when he was forty years old. He communicated this in a letter to Goldbach 
dated 6 May 1747. The proof relies on infinite descent, and proceeds in five steps; 
we state the first step from that proof below as we will be using it in the next 
paragraph: 

• The product of two numbers, each of which is a sum of two squares, is itself 
a sum of two squares. 

Given an OPN iV = p k m 2 , since p = 1 (mod 4), by Fermat's theorem we 
can write p as a sum of two squares. By Euler's first step above, p k can likewise 
be expressed as a sum of two squares, p k being the product of k p's. Hence, we 
can write p k = r 2 + s 2 for some positive integers r and s. Multiplying both sides 
of the last equation by m 2 , we get N = p k m 2 = m 2 (r 2 + s 2 ) = (mr) 2 + (ms) 2 . 
Hence, an odd perfect number may be expressed as a sum of two squares. 

Remark 3.4.1. Let 8(n) be the number of integers k < n that can be expressed 

6(n) 

as k = a 2 + b 2 , where a and b are integers. Does the limit lim exist and what 

is its value? Numerical computations suggest that it exists. 0(n) is approximately 

3 \ TXi 

= , so that the limit exists and equals zero. More precisely, the 



4/ \/hg(n) 

number of integers less than n that are sums of two squares behaves like K ■ 

Tl 

■. where K is the Landau- Ramanuj an constant. Dave Hare from the Maple 



\A>g(n) 

group has computed 10000 digits of K. The first digits are K w 0.764223653. 



See |http://www .mathsoft.com/asolve/constant/lr /lr.html| for more information. 
Therefore the sums of two squares have density 0. Since OPNs are expressible as 
sums of two squares, then OPNs have density 0, too. 



In the next chapter, we shall take a closer look into the nature of abundancy 
outlaws (which were first described in Section 2.2). We shall also describe a 
systematic procedure on how to bound the prime factors of an OPN N, using 
the latest current knowledge on N as well as some novel results. Lastly, we shall 
discuss some of the original results of the author pertaining to inequalities between 
the components of an OPN. 



Chapter 4 

OPN Solution Attempts 2: 
Some New Approaches 

In Chapter 3, we saw how increasing the lower bounds for ou(N) (the number 
of distinct prime factors of an OPN N) and N itself could potentially prove or 
disprove the OPN Conjecture. We also saw how the concept of divisibility may 
be used to derive congruence conditions for N. 

Here, we shall take a closer look into the following new approaches for at- 
tempting to solve the OPN Problem: 

• What are abundancy outlaws? How are they related to abundancy indices? 
How could one use the concept of abundancy outlaws to (potentially) dis- 
prove the OPN Conjecture? 

• How can one bound the prime factors of an OPN N? Is there a system- 
atic procedure on how to do this? (We discuss the author's results on the 
relationships between the components of N in Subsection 4.2.4.) 
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• Can we use the abundancy index concept to "count" the number of OPNs? 
(We answer this question in the negative for a particular case.) 

The reader is advised to review Section 2.2 of this thesis prior to commencing 
a study of this chapter. 



4.1 Abundancy Outlaws and Related Concepts 

Modern treatments of problems involving the abundancy index have been con- 
cerned with two fundamental questions: 



I. Given a rational number — , does there exist some positive integer x such that 

I(x) = ^ = %1 
x b 

II. When does the equation I(x) = — have exactly one solution for x? 

We give various answers to these two questions in the three subsections that 
follow. 



4.1.1 Friendly and Solitary Numbers 

If x is the unique solution of I(x) = — (for a given rational number — ) 
then x is called a solitary number. On the other hand, if x is one of at least two 
solutions of I(x) = — (for a given rational number — ) then x is called a friendly 
number. We formalize these two concepts in the following definition: 

Definition 4.1.1. Let x and y be distinct positive integers. If x and y satisfy 
the equation I(x) = I(y) then (x,y) is called a friendly pair. Each member of 



the pair is called a friendly number. A number which is not friendly is called a 
solitary number. 

We illustrate these concepts with several examples. 

Example 4.1.1. Clearly, if a and b are perfect numbers with a ^ b (i.e. cr(a) = 2a, 
a{b) = 2b), then (a, b) is a friendly pair. 

Example 4.1.2. We claim that, given a positive integer n satisfying gcd(n, 42) = 
1, (6n, 28n) is a friendly pair. To prove this, note that gcd(n,42) = 1 means 
that gcd(n, 2) = gcd(n, 3) = gcd(n, 7) = 1. Let us now compute J(6n) and 
J(28n) separately. Since gcd(n, 2) = gcd(n, 3) = 1, then gcd(n, 6) = 1 and 
J(6n) = J(6) J(n) = 2I(n) since 6 is a perfect number. Similarly, since gcd(n, 2) = 
gcd(n, 7) = 1, then gcd(n, 28) = 1 and 7(28n) = 7(28) I{n) = 2I(n) since 28 is a 
perfect number. These computations show that I(6n) = I(28n) = 2I(n) whenever 
gcd(n, 42) = 1, and therefore (6n, 28n) is a friendly pair for such n. 

Remark 4.1.1. Since there exist infinitely many positive integers n satisfying 
gcd(n, 42) = 1, Example 14.1.21 shows that there exist infinitely many friendly 
numbers. 

Example 4.1.3. M. G. Greening showed in 1977 that numbers n such that 
gcd(n,o~(n)) = 1, are solitary. For example, the numbers 1 through 5 are all 
solitary by virtue of Greening's criterion. There are 53 numbers less than 100, 
which are known to be solitary, but there are some numbers, such as 10, 14, 15, and 
20 for which we cannot decide "solitude" . (This is because it is, in general, difficult 
to determine whether a particular number is solitary, since the only tool that we 



have so far to make such determination, namely Greening's result, is sufficient 
but not necessary. In the other direction, if any numbers up to 372 (other than 
those listed in the Online Encyclopedia of Integer Sequences) are friendly, then 
the smallest corresponding values of the friendly pairs are > 10 30 [3T].) Also, we 
remark that there exist numbers such as n = 18,45,48 and 52 which are solitary 
but for which gcd(n, <r(n)) ^ 1. 

We give here a proof of Greening's criterion. Suppose that a number n with 
gcd(n, cr(n)) = 1 is not solitary. Then n is friendly, i.e. there exists some number 
x/n such that J(n) = I(x) . This is equivalent to xa(n) = na(x), which implies 
that, since gcd(n, cr(n)) = 1, n | x or x is a multiple of n. Thus any friend of n 
must be a (nontrivial) multiple of it (since n/i). Hence we can write x = mn 
where m > 2. Write m = jk with gcd(j, n) = 1 and gcd(k,n) > 1. Then by 
virtue of Lemma 12.2.21 I(x) > I(kn) since kn is a factor of x (unless j = 1), 
which implies that I(x) > I(n) (since k > 1 and this follows from gcd(fc, n) > 1). 
This is a contradiction. If j = 1, then we have x = kn with gcd(/c, n) > 1. Again, 
by virtue of Lemma I2.2.2l and similar considerations as before, I(x) > I(n) (unless 
k = 1, but this cannot happen since jk = m > 2) which is again a contradiction. 
Thus, numbers n with gcd(n, <r(n)) = 1 are solitary. 

Example 4.1.4. We claim that primes and powers of primes are solitary. It 
suffices to show that p k and o~{p k ) are relatively prime. To this end, consider the 
equation (p— l)a(p k ) = p k+l — l. Since this can be rewritten as (1— p)o~(p k )+p-p k = 
1, then we have gcd(p k ,a(p k )) = 1. By Greening's criterion, p k is solitary. Thus, 
primes and powers of primes are solitary. 



Remark 4.1.2. Since there are infinitely many primes (first proved by Euclid), 
and therefore infinitely many prime powers, Example 14.1.41 shows that there are 
infinitely many solitary numbers. 

Remark 4.1.3. While not much is known about the nature of solitary numbers, 
we do know that the density of friendly numbers is positive, first shown by Erdds 
[TT] : The number of solutions of 1(a) = 1(b) satisfying a < b < x equals Cx + o(x), 
where C > is a constant (in fact, C > j|= [2J). In 1996, Carl Pomerance 
told Dean Hickerson that he could prove that the solitary numbers have positive 
density, disproving a conjecture by Anderson and Hickerson in 1977. However, 
this proof seems not to ever have been published. 



4.1.2 Abundancy Indices and Outlaws 

On the other hand, rational numbers — for which I(x) = — has no solution 

b b 

for x are called abundancy outlaws. (Recall Definition 12.2.31 ) Of course, those 
rationals — for which I(x) = — has at least one solution for x are called abundancy 
indices. 



It is best to illustrate with some examples. 



Example 4.1.5. At once, Lemma 12.2.51 reveals a class of abundancy outlaws. 

n- i i . m . , . , m a(n) . . . 

Since that lemma says that — is an outlaw when 1 < — < (with gcd(m, n) = 

n n n 

1), then we have the class ^ — of outlaws (with t > 1). (We shall show later 

that, under certain conditions, - is also an abundancy outlaw.) 



Example 4.1.6. Let a,b,c be positive integers, and let p be a prime such that 

gcd(a,p) = 1, b = p c (so that b is a prime power, and gcd(a, b) = 1), and a > cx(6). 

<r(n) a 

Suppose we want to find a positive integer, n, such that = -. (That is, we 

no 

want to determine if — is an abundancy index or not.) This problem is equivalent 
to the problem of finding positive integers to, k such that: 

• n = mp k , k > c, and gcd(m,p) = 1 (or equivalently, gcd(m, b) = 1) 
a{m) ap k ~ c 



to a(p k ) 

We will formally state this result as a lemma later (where we will then present a 

proof), but for now let us see how we may apply this result towards showing that 
7 

the fraction - is an abundancy index. (Indeed, we are then able to construct an 

V " + 4.- C ■ a ( ■ a \ 

explicit n satisfying = - for a given -.) 

n b b 

ct(ti) 1 

We now attempt to find a positive integer n such that = -. This 

n 2 

problem, is equivalent to the problem of finding positive integers to, k such that: 

1.1 n = 2 k m, where gcd(2,TO) = 1 (that is, k is the largest power of 2 to divide 
<t(to) _ 7 • 2 k ~ 1 
to a(2 k ) 

We will now check different values of k, attempting each time to find to 
satisfying these conditions subject to the choice of k. For each k, we will proceed 
until one of the following happens: 

• We find to satisfying (1.1) and (1.2). In this case, n = 2 k m is a solution to 
our problem. 
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• We prove that there is no m satisfying (1.1) and (1.2). In this case, there is 
no solution to our problem of the form n = 2 k m, where gcd(2, m) = 1. 

• The problem becomes impractical to pursue. Often a given value of k will 
leave us with a problem which either cannot be solved with this method, or 
which is too complicated to be solved in a reasonable amount of time. 

We start with k — 1. Then our conditions are: 

2.1 n = 2m, and gcd(m, 2) = 1 
22 a(m) _ 7 • 2 1 - 1 _ 7 
m a(2 1 ) 3' 

Thus, our problem is to find m satisfying (2.1) and (2.2). Let us carry the 
process one step further for the case k — 1. To do this, we will treat m in the 
same manner in which we initially treated n. Our goal is to find mi, k\ such that: 

m = 3 kl m 1 , and gcd(mi,3) = 1. (Note that since m 1 \ m, and gcd(m, 2) = 1, 

we actually need gcd(mi,6) = 1.) 
o-(mi) _ 7-3 fel - 1 
mi a(3 kl ) 

Let us now check the case of k — 1, k\ — 1. Our goal is to find m x such 
that: 

4.1 m — 3mi, and gcd(mi, 6) = 1 

42 = 7 

mi 4 



3.1 
3.2 
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Thus, we want to find some positive integer m\ such that gcd(mi,6) = 1 

(j(tyi ~) 7 o~(m \ 7 

and = -. However, if = -, then 4 I m x (since gcd(4, 7) = 1). If 

mi 4 mi 4 

4 | mi, then gcd(mi,6) ^ 1, a contradiction. Therefore, there is no such mi. 



Our method has shown that, in the case of k — 1, k\ — 1, there is no positive 
integer n which solves our original problem. In particular, n is not of the form 
n = 2 1 m = 2 1 (3 1 (mi)) = 6mi, where gcd(mi, 6) = 1. (Another way to say this is 
that 6 is not a unitary divisor of any solution to our problem.) 



(Note that this does not prove the nonexistence of a solution to our problem; 
it only disproves the existence of a solution of the form given in the last paragraph. 
In order to disprove the existence of a solution of any given problem, we have to 
show that no solution exists for any value of k. Here we have not even eliminated 
the case of k — 1, but only the special case where k\ = 1.) 



Let us move on to k = 2. Our goal is to find m such that: 



5.1 n = 2 2 m = 4m, with gcd(m, 2) = 1 



q(m)_ 7-2 2 - 1 _14_ o 
b - Z m ~ <x(2 2 ) _ 7 



o~{ m ) 

Here, we must find m such that m is odd and = 2; that is, m must be 

m 

an odd perfect number. If m is an odd perfect number, then n = 4m is a solution 
to our problem. This is not especially helpful in our search for a solution, so we 
will move on to another case. 
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Here we will skip the cases k — 3 and k — 4, because they are not especially 
interesting compared to the next case we will deal with. 

Consider k — 5. Our goal now is to find m such that: 

6.1 n — 2 5 m = 32m, and gcd(m, 2) = 1 
g2 ff(m) 7-2 4 112 16 



m (x(2 5 ) 63 9 ' 

We can now apply our method to m. Keep in mind that the process will be 

16 

slightly different this time, since the denominator of — is a prime power, not just 

9 

a prime. Here, b = 9 = 3 2 , so we will use p = 3 and c = 2 (as they are used in the 
beginning of this example). Our goal is to find positive integers m l5 k x such that: 



7.1 m = mi • 3 fel , k± > 2, and gcd(mi,3) = 1. (Note that m x | m, and 

gcd(m, 2) = 1, so gcd(mi, 6) = 1.) 
72 tr(mi) _ 16 • 3 fcl ~ 2 
mi (j(3 fel ) 

We will consider two of the possible cases here: ki — 2 and ki = 3. 

First, let h = 2. We get 

c{mi) 16 16 
mi ~ (j(3 2 ) ~ 13' 

Carrying the process one step further, we will search for such an mi. We must 
find positive integers m 2 , fc 2 such that: 
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8.1 mi = m 2 • I3 k2 , and gcd(m 2 , 13) = 1. (Note that m 2 \ mi, 



so 



8.2 



gcd(m 2 , 6) = gcd(m 2 , 13) = 1.) 
a{m 2 ) _ 16 • 13* 2 " 1 
m 2 ~ cr(13 fc 2) 

Let k 2 = 1. Then mi = 13m 2 , and 

a(m 2 ) _ 16 _ 8 
m 2 ~~ 14 ~~ 7' 

Let m 2 = 7. Then m 2 satisfies both (8.1) and (8.2) of the case 
2, k 2 = 1. This gives us a solution; all we have to do now is work backwards until 
we get n. 

First, mi = 13m 2 = 13-7. Next, m = 3 kl m x = 13 • 7 • 3 2 . Finally, 
n = 2 5 m = 2 5 • 3 2 • 7 • 13 = 26208. 



(jf??/! 7 

Thus, we have found a solution to the problem = -. 

n 2 



Now we will consider the case k± = 3, which will give us one more solution: 

cx(mi) _ 16 • 3 _ 48 _ 6 
mi ~ (r(3 3 ) ~ 40 ~ 5' 

(j(tyi ") 6 

That is, we need to find mi such that = -, and gcd(mi,6) = 1. If we 

mi 5 

let mi = 5, then we have solved this problem, and have thus discovered another 

ct(ti) 7 

solution of the problem = -, this time for the case k — 5, k± — 3. Again, we 

n 2 

work backwards until we get n. 



First, m = 3 fcl mi = 3 3 • 5. Next, n = 2 5 m = 2 5 • 3 3 • 5 = 4320. 



This is a second solution to the problem 
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a(n) 7 



n 2 



7 

Remark 4.1.4. Example 14.1.61 shows that - is an abundancy index. Also, since 

7 

the equation I(n) = - has at least two solutions, namely rii = 4320 and n 2 = 
26208, this implies that rii and n 2 here are friendly. 



Prior to discussing the proof of the lemma outlined in Example I4.1.6[ we 
review some known properties of the abundancy index: 

Lemma 4.1.1. Properties of the Abundancy Index 

k 

• If a = Y\_Pj n3 ' w ^ ere PiiP2,P3, ■ ■ ■ ,Pk are distinct primes, k is a positive 

3=1 

integer, and the integral exponents n 1; n 2 , n 3 , . . . , n k are nonnegative, then 
1(a) = I I — —. -, where I is multiplicative. 

p 

• Ifp is prime then the least upper bound for the sequence {I{p n )}^ =0 is -. 

• If a% | a and ai > 0, then 1(a) > I(a\). 

• Obviously, if a\ \ a and 1(a) = I(ai), then a = a±. 

• If gcd(a,a(a)) = 1, then the unique solution to I(x) = 1(a) is x = a. 

Proof. Only the last assertion is not so obvious. Note that I(x) = 1(a) is equiv- 
alent to a ■ a(x) = x ■ a (a). If a and a(a) are coprime, then a | x (so that a 
divides every solution) . By the fourth result of Lemma 14.1.11 x = a is the sole 
solution. □ 
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We discuss some known properties that can help us decide whether a par- 

r 

ticular fraction - is an abundancy index: 
s 

r 

Lemma 4.1.2. When is - an abundancy index? 

s 

r 

For which rational numbers - will 

s 

i( X ) = - n 

s 

r 

have at least one solution? In order for (*) to have solutions, - must be greater 

s 

r 



than or equal to one. If - — I, then by Lemma 4-1-1. x — 1 is the unique solution. 

s 

Throughout the rest of this lemma, it will be assumed that r and s represent given 
positive integers which are relatively prime, and that r > s. Let us now state some 
known results: 

• Note that s must divide (every solution for) x in (*). This property is easy 

r 

to observe since I(x) = - implies that s ■ c(x) = r ■ x, and gcd(r, s) = 1. 

s 

r o~{ s ) 

• If a solution to (*) exists, then r > a(s) (since - = I(x) > I(s) = ). 

s s 

r 

• If 1(a) > -, then (*) has no solution which is divisible by a. Additionally, 

s 

if s is divisible by a, then (*) has no solution. 

• {1(b) : b G Z + } is dense in the interval (1, oo). 

T 

• The set of values - for which (*) has no solution is also dense in (1, oo). 

Proof. Here, we prove the third assertion. Suppose (*) has a solution which 

is divisible by a G Z + . Then a \ x, and by the third result in Lemma 14.1.11 

I(x) > 1(a). But then, by assumption 1(a) > -, which implies that 1(a) > - = 

s ' s 

I(x) > 1(a), resulting in the contradiction 1(a) > 1(a). Thus, (*) has no solution 
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which is divisible by a. The second statement follows from this and the first result 
in Lemma [4.1.21 (i.e. s must divide every solution for x). □ 

We now state and prove the following lemma (taken from [17]) which was 
used in Example 14.1.61 

Lemma 4.1.3. (Ludwick) Let a,b,c G Z + ; and letp be a prime such that gcd(a,p) = 
1, b = p c (so that b is a prime power, and gcd(a, b) = 1), and a > o~(b). Sup- 
pose we want to find a positive integer, n, such that - - = y. This problem is 

n b 

equivalent to the problem of finding positive integers m, k such that: 

• n = mp h , k > c, and gcd(m,p) = 1 (or equivalently, gcd(m, b) — 1) 
a{m) ap k ~ c 



m o~(p k ) 

Proof. First, we will show that if we can find n G Z + satisfying - - = — , then 

n b 

we can find m,k G Z + satisfying the two itemized conditions above. 



Suppose we have n G Z + such that — = — . Then by Lemma [4. 1.2\ b \ n; 

no 

that is, p c | n. Since 6 is a prime power, there is some k G Z + such that n = mp k , 
with gcd(m, 6) = 1. Here, k is the largest power of p that divides n, so clearly 
k > c. This satisfies the first condition. Now, we have 

a{n) a a 
n bp 

and we also have 



o~(n) a(p k ) a{m) 

Thus, 



n p k m 
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a a{p k ) a(m) 

and so 



pc pk m 



a(m) ap 



k—c 



m 



a{p k ) 



This satisfies the second condition. Therefore, solving = y for n gives us 

n b 

m,k G Z + satisfying the two conditions. 



Conversely, we will show that if we can find m, k G Z + satisfying the two 

conditions, then we can find n G Z + such that = — . 

n b 

Suppose we have m,k G Z + satisfying the two conditions. Let n = mp k . 

Then, 

a(n) a(m)a(p k ) ap k ~ c a(p k ) ap k ~ c a a 



n m p k a{p k ) p k p k p c b 



□ 



4.1.3 OPNs, Abundancy Outlaws and the Fraction 



p + 2 



P 

After defining the abundancy index and exploring various known properties, 
we briefly discuss some related concepts. Positive integers having integer-valued 
abundancy indices are said to be multiperfect numbers. One is the only odd 
multiperfect that has been discovered. Richard Ryan hopes that his "study of 
the abundancy index will lead to the discovery of other odd multiperfects" , or to 
the proof of their nonexistence. Since the abundancy index of a number n can 
be thought of as a measure of its perfection (i. if I(n) < 2 then n is deficient; 
ii. if I(n) = 2 then n is perfect; and iii. if I(n) > 2 then n is abundant), it is 
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fitting to consider it a very useful tool in gaining a better understanding of perfect 
numbers. In fact, Judy Holdener [32] proved the following theorem which provides 
conditions equivalent to the existence of an OPN: 

Theorem 4.1.1. There exists an odd perfect number if and only if there exist 
positive integers p, n and a such that p = a = 1 (mod 4), where p is a prime not 

A A' A Tf \ 2 P >-1) 

dividing n, and l[n) = — — . 

Proof. By Euler's characterization of an OPN N = p a m 2 , it must be true that 
p is a prime satisfying gcd(p, m) = 1 and p = a = 1 (mod 4). Hence c(N) = 
a(p a m 2 ) = a(p a )a(m 2 ) = 2p a m 2 , and 

T , 9N aim 2 ) 2p a 2p a (p-l) 

1 [771 ) = — — = = 

This proves the forward direction of the theorem. 

„ , ... , , r/ \ 2p°(p — 1) 

Conversely, assume there is a positive integer n such that 1 (n) = — — , 

where p = a = 1 (mod 4) and p is a prime with p \ n. Then 

7(n • p a ) = I(n) ■ I(p a ) = 2 P a [P~y P °^~\ = 2. 
v ' K ' K ' p a+1 - 1 p a {p- 1) 

So n ■ p a is a perfect number. 

Next, we claim that n • p a cannot be even. Suppose to the contrary that 
n ■ p a is even. Then it would have the Euclid-Euler form for even perfect numbers: 

n-p a = 2 m " 1 (2 m - 1) 

where 2 m — 1 is prime. Since 2 m — 1 is the only odd prime factor on the RHS, 
p a = p 1 = 2 m — 1. But p = 1 (mod 4) and 2 m — 1 = 3 (mod 4) (because m must 
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be at least 2 in order for 2 m — 1 to be prime). This is clearly a contradiction, and 
thus n ■ p a is not even. Consequently, n ■ p a is an OPN. □ 

By Theorem I4.1.1[ it follows that if one could find an integer n having abun- 
5 

dancy index equal to - (which occurs as a special case of the theorem, specifically 
for p = 5 and a = 1), then one would be able to produce an odd perfect number. 
Here we then realize the usefulness of characterizing fractions in (l,oo) that are 
abundancy outlaws. (Recall Definition 12.2.31 ) 

Let us now consider the sequence of rational numbers in (l,oo). (Note 
that, since the number 1 is solitary and 1(1) = 1, the equation I(x) = 1 has the 
lone solution x = 1.) For each numerator a > 1, we list the fractions — , with 
gcd(a, b) = 1, so that the denominators 1 < b < a appear in increasing order: 

23344555566777777 



1' 1'2' 1'3' 1'2'3'4' 1'5' 1'2'3'4'5'6'"' 
It is intuitive that each term in this sequence must be either an abundancy index 
or an abundancy outlaw, but it is, in general, difficult to determine the status of 
a given fraction. We may thus partition the sequence into three (3) categories: 
(/) those fractions that are known to be abundancy indices, (//) those that are 
known to be abundancy outlaws, and (///) those whose abundancy index/outlaw 
status is unknown. We wish to capture outlaws from the third category, thereby 

increasing the size of the second category. Since fractions of the form - — 

for t > 1 belong to the first category (by Lemma 12.2.51 and Example 14.1. 5ft . 
it is tempting to consider fractions of the form ~~^y~~- Judy Holdener and 
William Stanton proved in 2007 [51] that, under certain conditions, ^ - is 
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an abundancy outlaw. They noted that their original interest in such fractions 

stemmed from their interest in the fraction - = - - — . Unfortunately, the 

3 3 J ' 

results they obtained do not allow them to say anything about fractions of the 

form a ^ — — = — . Such elusive fractions remain in category three. 
p p 

Equivalently, we may ask: Does there exist an odd number s G Z + (with 

s + 2 

s > 1) such that I(x) = has at least one solution? The answer to this 

s 

question is unknown, but we can go ahead and discuss some properties. 



By Lemma f4.1.2[ s \ x. Assume s is an odd composite; then a(s) > 1 + s + d, 

where d is a divisor of s satisfying 1 < d < s. Since s is odd, d > 3, which means 

s + 2 

that ens) > s + 4 > s + 2, or I(s) > , a contradiction. Hence, s must be 

s 

prime. If 1 < c < s, then gcd(c, s) = 1 (since s is prime), and we have: 

I(cs) = 1(c) ■ I s > > ^ ^- = + - > . 

c s s z s z \s J s 

Thus, x does not have a factor between 1 and s. Moreover, x is a perfect square; 
otherwise cr(x) would have a factor of 2 that cannot be "canceled" since the 
denominator, x, is odd. (For the same reason, whenever r and s are both odd, 

T 

any odd solution to I(x) = - must be a perfect square.) We also claim that 

s 

s + 2 

(s + 2) \ x and we prove this by showing that / {s 2 (s + 2) 2 ) > : 



, 9 , .,v 1 1 1 1 1 1 

/(s 2 (s + 2) 2 ) > 1+-+^+ -+ 7 -g+ — +- 



s s 2 s + 2 (s + 2) 2 s(s + 2) s 2 {s + 2) s{s + 2) 2 s 2 (s + 2) 2 



s 2 + s + 1 1 
7 2 + ^+2 



1 1 
1 + + - 

s + 2 s 



1 

+ 



s(s + 2) 



1 1 1 

s + s + 2 + s(s + 2) 



s 2 + s + 1 , s 2 + 4s + 2 2s + 3 



s 2 s(s + 2) 2 s 2 (s + 2) 2 
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(s 2 + s + 1) (s + 2) 2 + s(s + 2) 2 + 3 
s 2 (s + 2) 2 

(s + l) 2 (s + 2) 2 + 3 
s 2 (s + 2) 2 

s + l\ 2 3 s 2 + 2s + l s + 2 

s / + s 2 (s + 2) 2 > s 2 > s 



Using some of the principles in the last paragraph, Richard Ryan wrote a 

s + 2 

simple computer program which verified that I(x) = has no solution less 

s 

,16 



than 10 (when s G Z + is odd with s > 1) [58] . 



We now state (without proof) the conditions obtained by Holdener and 
Stanton [3l] in order for the fraction - to be an abundancy outlaw. 

Theorem 4.1.2. For a positive integer t, let - be a fraction in lowest 



terms, and let N = ~Y\_Pi kl f or V r ^ mes P11P21 ■ ■ ■ iPn- If there exists a positive 

i=i 

integer j < n such that pj < \o~{-pq) and a {Pj ki ) has a divisor D > 1 such that 



at least one of the following is true: 

1. I{p^)I{D) > ff(A ^ +t and gcd(At) 

2. gcd(£>, Nt) = 1 

, o-(N)+t , , 

then — is an abundancy outlaw. 



The following are some sequences of abundancy outlaws which can be con- 
structed from Theorem 14.1.21 
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• For all natural numbers m and nonnegative integers n, and for all odd primes 

a(2 m p 2n+1 ) + 1 

p such that gcd(p, cr(2 m )) = 1, the fraction — 2n+1 is an abundancy 

2 m p n 

outlaw. 

• For all primes p > 3, - ^ — is an abundancy outlaw. If p = 2 or p = 3 

2p 

then *w + 1 is an ab Mex 

2p 

• If iV is an even perfect number, - 2p^~ — * s an a bundancy outlaw. 

• Let M be an odd natural number, and let p, a, and t be odd natural numbers 

a(p a M) + 1 

such that p \ M and p < \a(M). Then, if is in lowest terms, 

* p a M 

g ( p ° A 2 + ' » *-d-y outlaw. 

• For primes p and q, with 3 < q, p < q, and gcd(p, q + 2) = gcd(g, p + 2) = 1, 

^ g) + 1 is an abundancy outlaw. 
pq 

Remark 4.1.5. The last assertion in the preceding paragraph produces outlaws 
with ease. We illustrate this using odd primes p and q satisfying 3 < p < q and 
q = 1 (mod p). It follows that p \ (q + 2) and q \ (p + 2). By Dirichlet's theorem 
on arithmetic progressions of primes, we are assured of the existence of an infinite 
sequence of primes q satisfying q = 1 (mod p). Thus, there is an infinite class of 
abundancy outlaws corresponding to each odd prime p > 3. 



Judy Holdener, with Laura Czarnecki, also obtained the following results in 
the summer of 2007 [35] : 
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Theorem 4.1.3. If — is a fraction greater than 1 in reduced form, — = I(N) 

b b 



for some N G N, and b has a divisor D = J^Pi fe4 such that I(piD) > — for all 

i=l 



1 < i < n, then — — -— is an abundancy index as well. 
a(D) b 

Proof. Suppose that I(N) = — for some AT e N. Then by Lemma I4.1.2[ b I A" 

b 

since gcd(a, b) = 1. Because I{piD) > — for all 1 < i < n, we know that it is 
impossible that piD | A" by the contrapositive of the third result in Lemma 14.1.11 
However, we know (by Lemma I4.1.2P that D | N, so pi ki+1 \ N . Thus we may 

/ n \ 

write = | J • r = Dr, where gcd(pj,r) = 1 for all 1 < i < n, that is, 

a a(N) 



gcd(Z^,r) = 1. Then, since o~(N) is multiplicative, we may write 



b N 

— — - — . Therefore, I Or) = - - = — r. Thus, if I(N) = — for some 
D r r b ct{D) b 

iVeN, then I Or) = —. — - for some r G N. □ 
b o\D) 

o~ (mn) -\- o~ (m)t 

Corollary 4.1.1. Let m,n,t G N. If is in reduced form with 

mn 

nk it/ cr(mn) + aOm)t a(n)+t 
Pi ' and Ivpim) > for all 1 <i < I, then is an 

mn ' n 

i=i 

abundancy index if — °^ — — is an index. 

mn 

Proof. The proof is very similar to that of Theorem l4.1.31 We only need to observe 
that, under the assumptions given in this corollary, gcd(m, n) = 1. □ 

Remark 4.1.6. If, in Corollary I4.1.1[ we have t — 1 and n = p for some prime 

p, then the corollary tells us that it is an abundancy index, then 

mp 

0{P) + 1 m , r 27 39 45 , 57 n i, 1 • r 

is as well. The tractions — , — , — and — all illustrate this fact. It we 

p 14' 22' 26 34 

9 13 15 

could determine that these are indeed indices, then we could say that -, — , — 

7 11 13 
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and — , all of the form — — , are indices as well. Fractions of the form ' y 



17 p p 

continue to elude characterization as indices or outlaws! This is significant because 

5 

Paul Weiner [73] proved that if there exists an integer N with abundancy — , then 

5iV is an odd perfect number. The fraction - is of the form — — for p — 3. 

3 p 

Theorem 4.1.4. If p > q > 2 are primes satisfying p > q 2 — q — 1, then 
cr(qp) +5 — I 



qp 



is an abundancy outlaw. 



o~(qp) + q — 1 

Proof. Suppose that I(N) = for some N G N. Because a{qp) + q — 

qp 

1 = (g + l)(p + 1) + q — 1 = qp + 2q + p and p > q > 2 are primes, we know that 

gcd(a(qp) + q — 1, qp) = 1, and by Lemma [4.1.21 we have qp \ N. If q 2 p | N, then 
a(qp) + q - 1 (r(q 2 p) 

> — - — . bimphtymg, this becomes 

qp q 2 p 

(g + i)(p+ i) + q- 1 > (g 2 + Q + 1)(P+ 1) 

and then g 2 — g — 1 > p. This contradicts our hypothesis that p > q 2 — q — 1, and 

therefore g 2 p f iV. Consequently, we may write N = qK, where gcd(g, K) = 1. 

Then, we have I(N) = Uq)I(K), which gives UK) = = G ^ qP) + q ~ l . 

7(g) gp 

g+1 

g p ' ~~ o ^ g <?+! 9 + 1 

= , 1 — . On setting m = p \-q and n = p and observing 

g+1 g+1 6F 2 y F 2 6 

that n < m < o~(n) with gcd(m, n) = 1, then by Lemma f2 . 2 . 51 ( and Example l4.1.5p . 

m c(qp) + g — 1 g . 

— = • is an abundancy outlaw. But this contradicts the tact 

n qp g+1 

I(N) <j(qp) + g — 1 g 

that I(K) = — -— = ■ for some K E N, whence it follows 

/(g) qp g+1 

that I(N) ^ a ^ + q - 1 for all JVeN. Consequently, + 9 ~ 1 must be 

qp qp 

an abundancy outlaw under the conditions specified for the primes p and g. □ 



These results allow us to move a few more fractions in (1, oo) from the set 

of infinitely many fractions that we are unable to classify (category III), into the 

infinite set of fractions that are abundancy outlaws (category II). Furthermore, we 

can see that certain fractions are linked to others in important ways: determining 

the status of a given fraction can lead to the classification of new abundancy 

outlaws and indices. If the converse of Theorem 14.1.31 could be proved, then we 

would be able to divide certain fractions into equivalence classes of sorts, that 

is, sets of fractions with the same abundancy index/outlaw status. However, the 

question of the existence of an OPN (e.g. the status of the fraction — — for 

P 

an odd prime p) remains as elusive as ever. 



4.2 Bounds for the Prime Factors of OPNs 



In this section, bounds for each of the distinct prime factors of an OPN N are 
derived, drawing heavily from existing works. We do this using cases based on the 
total number of distinct prime factors of N (i.e. fl(iV)). We also study further 
results in the field and give examples of various techniques used, including an in- 
depth and detailed discussion of the factor chain approach. We give in Subsection 
4.2.3 explicit double-sided bounds for each of the prime factors of an OPN N with 
co(N) = 9. We end the section with a discussion of the author's results on the 
relationships between the components of an OPN N. 
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4.2.1 Results on OPNs 

Let iV = qi ai q2 a2 • • -qt at be the canonical factorization of an OPN N (i.e. 
qi, q2, ■ ■ ■ , qt are distinct primes with qi < q 2 < . . . < qt and t = u(N)). Then the 
following statements are true: 

• Qt > 100000007 from Goto and Ohno, improving on Jenkins 

• qt-i > 10007 from Iannucci, improving on Pomerance 

• Qt-2 > 101 from Iannucci 

• qi < 2 2% 1 (t — i + 1) for 2 < i < 6 from Kishore 

• q\ < from Grim 

• qk ak > 10 20 for some k from Cohen, improving on Muskat 

• N > 10 300 from Brent, et. al, improving on Brent and Cohen (A search is 



currently on in |http://www.oddperfect.org to prove that N > 10 500 . 



• N = 1 (mod 12), N = 81 (mod 324) or iV = 117 (mod 468) from Roberts, im- 
proving on Touchard and Holdener 

• N < 2 4 ' from Nielsen, improving on Cook 

t 

• = n(N) > 75 from Hare, improving on Iannucci and Sorli 

i=l 

• t > 9 from Nielsen, improving on Hagis and Chein 

• t > 12 if q\ > 5 from Nielsen, improving on both Hagis and Kishore 

• t > 17 if q\ > 7 from Greathouse, improving on Norton 

• t > 29 if qi > 11 from Greathouse, improving on Norton 



Suryanarayana and Hagis [27J showed that, in all cases, 0.596 < J2 p \n ~ < 
0.694. Their paper gives more precise bounds when A" is divisible by 3 or 5 (or 
both). Cohen [7J also gave more strict ranges for the same sum, including an 
argument that such bounds are unlikely to be improved upon significantly. 

We use the preceding facts about OPNs to derive explicit double-sided 
bounds for the prime factors of an OPN N with u(N) = 9, in Subsection 4.2.3. 

4.2.2 Algorithmic Implementation of Factor Chains 

(We borrow heavily the following material from [63].) 

In the discussion that follows, we will let iV denote an OPN, assuming one 
exists, with the prime decomposition 

U V w 

i=l i=l i=l 

which we interpret as follows: each pi ai is a known component of N, each is a 
known prime factor of N but the exponent fe, is unknown, and each prime factor 
r, of iV and exponent q are unknown. By "known", we mean either explicitly 
postulated or the consequence of such an assumption. Any of u, v, w may be 
zero, in which case we set A, (/,, u, respectively, equal to 1. We also let fn denote 
a proper divisor of m (except 1 = 1). 

We can now illustrate the factor chain approach via an algorithmic imple- 
mentation that can be used to test a given lower bound for u(N). We assume 
that is an OPN with u(N) = t distinct prime factors. In brief, the algorithm 



may be described as a progressive sieve, or "coin-sorter", in which the sieve gets 
finer and finer, so that eventually nothing is allowed through. We shall use the 
terminology of graph theory to describe the branching process. Since 3 | N if 
t < 11, for our present purposes the even powers of 3 are the roots of the trees. 
If 3 2 is an exact divisor of N, then, since cr(N) = 2N, a(3 2 ) = 13 is a divisor 
of N, and so the children of the root 3 2 are labelled with different powers of 13. 
The first of these is 13 1 , meaning that we assume that 13 is an exact divisor of N 
(and hence that 13 is the special prime), the second 13 2 , then 13 4 , 13 5 , . . .. Each 
of these possibilities leads to further factorizations and further subtrees. Having 
terminated all these, by methods to be described, we then assume that 3 4 is an 
exact divisor, beginning the second tree, and we continue in this manner. Only 
prime powers as allowed in Subsection 4.2.1 are considered, and notice is taken 
of whether the special prime has been specified earlier in any path. These powers 
are called Eulerian. 

We distinguish between initial components, which label the nodes and com- 
prise initial primes and initial exponents, and consequent primes, which arise 
within a tree through factorization. It is necessary to maintain a count of the 
total number of distinct initial and consequent primes as they arise within a path, 
and we let k be this number. 

Often, more than one new prime will arise from a single factorization. All 
are included in the count, within k, and, whenever further branching is required, 
the smallest available consequent prime is used as the new initial prime. This 
preferred strategy will give the greatest increase in I(Xp). On the other hand 



a strategy of selecting the largest available consequent prime will usually give a 
significant increase in k. 

To show that t > u>, say we build on earlier results which have presumably 
shown that t > u> — 1, and we suppose that t — u> — 1. (The reader may want 
to review Section 3.f at this point.) If, within any path, we have k > u> — 1, 
then there is clearly a contradiction, and that path is terminated. This is one of 
a number of possible contradictions that may arise and which terminate a path. 
The result will be proved when every path in every tree has been terminated 
with a contradiction (unless an OPN has been found). The different possible 
contradictions are indicated with upper case letters. 

In the contradiction just mentioned, we have too Many distinct prime factors 
of N: this is Contradiction Ml. If there are too Many occurrences of a single prime 
this is Contradiction M2; that is, within a path an initial prime has occurred as 
a consequent prime more times than the initial exponent. (So counts must also 
be maintained within each path of the occurrences of each initial prime as a 
consequent prime.) 

If k = u> — 3 but none of these k primes exceeds 100, then Iannucci's result 
must be (about to be) violated: this is Contradiction P3. If k = u — 2 and none of 
these primes exceeds 10 4 , then again, Iannucci's result is violated: Contradiction 
P2. Or if, in this case, one exceeds 10 4 but no other exceeds 100, then this is 
another version of Contradiction P3. If k — ou — 1 and none of these primes 
exceeds 10 8 , then Goto/Ohno's result is violated: Contradiction PI. In this case, 



there are the following further possibilities: one prime exceeds 10 but no other 
exceeds 10 4 , or one exceeds 10 8 , another exceeds 10 4 , but no other exceeds 100. 
These are other versions of contradictions P2 and P3, respectively. These, and 
some of the other forms of contradiction below, require only counts or comparisons, 
and no calculations. 

At the outset, a number B is chosen, then the number of subtrees with a 
given initial prime p is bounded by taking as initial components Eulerian powers p a 
with p a+1 < B. If possible, these trees are continued by factorizing a(p a ). When 
a becomes so large that p a+l > B, which may occur with a = 0, then we write q b 
for p a and we have one more subtree with this initial prime; it is distinguished by 
writing its initial component as q°°. This tree must be continued differently. In 
the first place, the smallest available consequent prime, which is not already an 
initial prime, is used to begin a new subtree. If no such primes are available, then 
we opt to use the procedure that follows. 

The product of the u initial components p a within a path is the number 
A. Those initial primes q with exponents oo, and all consequent primes which 
are not initial primes, are the v prime factors of /i. If k < ou — 1 then there are 
w = oj — k — 1 remaining prime factors of N, still to be found or postulated. 
These are the prime factors r of v. The numbers u, v, w are not fixed; they vary 
as the path develops, for example, by taking a consequent prime as another initial 
prime. 
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If factorization can no longer be used to provide further prime factors of 
N, so, in particular, there are no consequent primes which are not initial primes, 
then the following result (with proof omitted) is used: 

Lemma 4.2.1. Suppose w > 1, and assume r± < r-i < ■ ■ ■ < r w . Then 

I{XjLr x c ^ x ) 

< r 

2-I(Xp) ~ 

for r = r\, with strict inequality if v > 1 or w > 2. Further, if I^X^i 00 ) < 2, then 

2 + I(Xfj,°°)(w- 1) 
r < 2 - /(A/1 00 ) 

for r = ri . 

Here, Jjl is taken to be the product of powers q@, where q \ /i and (3 is given as 
follows. Let 60 — m i n {b '■ q h+l > -B}- If &o — 0, then we proceed in a manner to 
be described later. Otherwise, let 







bo, if bo is even (bo > 0), 

6 + l, if bo is odd, 



with one possible exception. If tt is the special prime, tc \ A and the set Q\ = 
{q : q = b = 1 (mod 4)} is nonempty, then take (3 = b for q = min Qi- Values of 
I(p a ) and I(qP) must be maintained, along with their product. This is the value 
of I(Xp) to be used in the result mentioned in this paragraph. We shall refer to 
Lemma 14.2.11 as Lemma X. 



Lemma X is used to provide an interval, the primes within which are con- 
sidered in turn as possible divisors of v. If there are No primes within the interval 



that have not been otherwise considered, then this is Contradiction N. New primes 
within the interval are taken in increasing order, giving still further factors of N 
either through factorization or through further applications of Lemma X. There 
will be occasions when no new primes arise through factorization, all being used 
earlier in the same path. Then again Lemma X is used to provide further possible 
prime factors of N (or, if k — u — 1, we may have found an OPN). This lemma 
specifically supplies the smallest possible candidate for the remaining primes; a 
still Smaller prime subsequently arising through factorization gives us Contradic- 
tion S. 

We also denote by q any consequent prime which is not an initial prime, 
and, for such primes, we let Q 2 = {q : q = 1 (mod 4)}. Then, for such primes, we 
let P — 2 with the possible exception that, considering all primes q, we let (3 = b 
or 1, as relevant, for q = min (Q 1 U Q 2 ), if this set is nonempty. Again, the value 
of /(A/2), defined as before, must be maintained. If this value exceeds 2, we have 
an Abundant divisor of N, and the path is terminated: Contradiction A. This 
may well occur with k < u — 1. Values of I(q°°) must also be maintained. These, 
multiplied with the values of I(p a ), give values of /(A/i°°). If this is less than 2 
and k = uj — 1 then, for all possible values of the exponents b, the postulated 
number N is Deficient: Contradiction D. 

Contradictions A and D are in fact contradictions of the following lemma, 
which we shall refer to as Lemma Y: 

Lemma 4.2.2. For any OPN N = Xfiu, as given as before, we have I(Xp) < 2 < 
I(\pfj,°°). Both inequalities are strict if v > 0; the left-hand inequality is strict if 



w > 0. 

If, on the other hand, we have a postulated set of prime powers p a and q b , 
for which /(A/2) < 2 < J(A/i°°), then the main inequality in Lemma Y is satisfied 
and we have candidates for an OPN. If v = w = 0, so that we are talking only 
of known powers p a , then their product is an OPN. Our sieving principle arises 
when v > 0. 

In every such case where we have a set of prime powers satisfying the main 
inequality of Lemma Y, with v > 0, we increase the value of B and investigate 
that set more closely. With the larger value of B, some prime powers shift from 
/i to A, and allow further factorization, often resulting quickly in Contradiction 
Ml or S. The value of I{jx) increases, so the interval given by Lemma X shortens, 
and hopefully the case which led to our increasing B is no longer exceptional, 
or Contradiction A or D may be enforced. In that case, we revert to the earlier 
value of B and continue from where we are. Alternatively, it may be necessary 
to increase B still further, and later perhaps further again. When w — 0, since 
I{p) — > as B — > oo, such cases must eventually be dispensed with, one way 

or the other. 



We summarize the various contradictions in the following table: 
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A 


There is an Abundant divisor. 


D 


The number is Deficient. 


Ml 


There are too Many prime factors. 


M2 


A single prime has occurred too Many times (an excess of that prime). 


N 


There is no New prime within the given interval. 


PI 


There is no Prime factor exceeding 10 8 . 


P2 


There is at most one Prime factor exceeding 10 4 . 


P3 


There are at most two Prime factors exceeding 100. 


S 


There is a prime Smaller than the purportedly smallest remaining prime. 


n 


None of the primes can be the special prime. 



One of these, Contradiction n, was not discussed previously. Within any 
path with k — u> — 1, if 7r is not implicit in an initial component and if there is no 
prime q = 1 (mod 4), then Contradiction IT may be invoked. 

Remark 4.2.1. Notice that the algorithm as presented in this subsection could 
be programmed directly to run on a high-speed computer (even desktop PCs). 
Prior experience with such algorithmic programs, however, has shown that it can 
take months (or years even), to check and/or test a particular value for t = u(N). 
Current computer architecture limits our capability to carry out these tasks at a 
reasonable amount of time. 

Example 4.2.1. Sigma chains (otherwise known as factor chains) are an easily 
automated system for proving facts about OPNs. Each line of the proof starts with 
a prime factor known or assumed to divide an OPN N, along with its exponent. 
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Since a is multiplicative, knowledge of this prime power leads to knowledge of 
other prime powers of N. If an impossibility arises (see below), that chain of the 
proof is terminated and the next possibility is considered. 

The following is the start of a proof that no OPN has a component less than 
10 30 . It would take many thousands of pages to complete this proof (and this 
has currently not been completed); this merely serves as an example of how 
one constructs such proofs. (The best-known result in this direction is that of 
Cohen [8]: An OPN has a component bigger than 10 20 .) 

The factor chains are terminated (the succeeding line is not to be indented 
further than the preceding) if it fails in one of the following ways: 

xs: The indicated prime appears more times than it is allowed, (e.g. If the 
chain assumes that 3 6 ||iV then a chain with 7 or more factors of 3 is terminated.) 

overabundant: The abundancy of the prime factors already exceeds 2, so 
regardless of the other factors, N will fail to be perfect. 

The factorizations of the largest half-dozen composites are due to the WIMS 
(WWWInteractive Multipurpose Server) Factoris at wims.unice.fr. (This proof is 
taken from [2TJ].) 



3° - > 1093 

1093 - > 2 * 547 

547 2 - > 3 * 163 * 613 

163 2 - >3*7*19*67 
7 2 - > 3 * 19 

19 2 - > 3 * 127 

127 2 - > 3 * 5419 

5419 2 - > 3 * 31 * 313 * 1009 
31 2 - > 3xs * 331 

31 4 — > 5 * 11 * 17351 overabundant 
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31 6 - > 917087137 

917087137 2 - > 841048817767943907 = 3xs * 43 * ... 

917087137 4 - > 707363113097541065394066657400343621 = 31747594185191 * 2228 . . . 9731 
31747594185191 2 - > 1007 . . . 1673 = 2671 * ... 

2671 2 - > 3xs * 7 * 19 * 31 * 577 

2671 4 - > 5 * ll 2 * 571 * 147389551 overabundant 

2671 6 - > 127 * 2860238405785894351 

2860238405785894351 2 - > 8180 . . . 1201 = 3xs * ... 
2860238405785894351 4 - > 
66928167 . . . 72862401 

= 5 * 11 * 27362961781 * ... overabundant 

2860238405785894351 6 - > 5475369 . . . 

. . . 5453601 = 7xs * 2339 * 337498477 * 1013 . . . 2827 * ... 

31747594185191 4 - > 1015882037027398808700619554107312555810842320401403361 

= 5 * 11 * 27581 * ... overabundant 

31747594185191 6 - > 

1023917 . . . 0882641 = 29 * 68279 * 17581747 * ... 

29 2 - > 13 * 67 

13 2 - > 3xs * 61 
13 4 - > 30941 



4.2.3 Explicit Double- Sided Bounds for the Prime Factors 

The results from Subsection 4.2.1 give some restrictions on the magnitude 
of the prime factors of an OPN N. For instance, we saw from Subsection 3.5.2 
that N is not divisible by 3 ■ 5 • 7 = 105. Consequently, it must be true that the 
third smallest prime factor g 3 > 11. To further derive bounds for the other prime 
factors, we will use some of the many published results on OPNs, a compendium 
of which has been presented in Subsection 4.2.1. 

For the largest prime factors of an OPN, Iannucci and Jenkins have worked 
to find lower bounds. The largest three factors must be at least 100000007, 10007, 
and 101. Goto and Ohno verified that the largest factor q t must be at least 
100000007 using an extension of the methods of Jenkins. 
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Nielsen, improving the bound of Hagis and Kishore, showed that if an OPN 
is not divisible by 3, it must have at least 12 distinct prime factors. Nielsen also 
showed that a general odd perfect number, if it exists, must have at least 9 distinct 
prime factors. Therefore, for 9 < t < 11, we have qi = 3. 

A result by Gran (and perhaps, independently too, by Perisastri) will be 
useful for our purposes later: q± < |t + 2. Results similar to those previously 
mentioned reduce the practicality of Gran's findings. In fact, a paper by Norton 
published about two years after supersedes Gran's inequality, except that Norton's 
method is slightly more computationally intensive. 

With an application of Goto/Ohno's and Iannucci's results, we can mod- 
estly improve on these bounds with an otherwise straightforward utilization of 
the abundancy index function: 

Let qi > 7. Now, suppose that t = u>(N) = 17. Using the lower bounds 
indicated for the three largest prime factors of N as before, we have: 

qi > 7 
q 2 > 11 

qs > 13 
q 4 > 17 

95 > 19 
q 6 > 23 
q 7 > 29 
qs > 31 
qg > 37 
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910 > 41 

911 > 43 

912 > 47 

913 > 53 

914 > 59 

915 > 101 
916 > 10007 

g 17 > 100000007 



Recall from Lemma [3.1.11 that 2 = — - — - < TT ( — — — ) . Also, note that 

Qi &i 

qi > di for all i implies that < — for all i, 

qi-l Oi — 1 



so that we have 



a ( N ) TT ( & \ / TT ( a 



=t ^ i=i 

The numerator of the rightmost fraction is approximately 

6.4778249375254265314282935191886 • 10 33 , while the denominator is approxi- 
mately 3.2172767985350308489460711424 • 10 33 , which gives a ratio of approxi- 
mately 

2.0134496790810999533771971435786, which is larger than 2. No contradiction at 
this point. 

As before, let qi > 7, but now suppose that t = u(N) = 16. Proceeding 
similarly as before, we have: 

91 > 7 

92 > 11 
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Again, we have 



9s > 13 

94 > 17 

95 > 19 

96 > 23 

97 > 29 

98 > 31 

99 > 37 

910 > 41 

911 > 43 

912 > 47 

913 > 53 

914 > 101 
915 > 10007 

916 > 100000007 



N 1A - \qi-lj f V a,; — 1 
where % > csj for each i. Our computations show that: 



1 ' 1.0979364300890553443098802574896 • 10 32 



n 



v Oi-l/ 55470289629914324981828812800000 
= 1.9793234133339626660318209208061 < 2. 

This results in the contradiction 2 < 2. We therefore conclude that t > 17 if gi > 7. 
We may likewise prove, using the same method, that t > 29 if q\ > 11. 



A result of great utility here is an earlier work of Kishore [44], where he 
proves that 
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qi < 2 T ~\t-i + 1) for2<i<6. 

These results (by Gran/Perisastri and Kishore) allow us to give explicitly reduced 

bounds for the lowest six (6) prime factors for an OPN N with a given number 

t = u(N) of distinct prime factors. For example, an OPN with nine (9) distinct 

2 • 9 

divisors has qi < 7 (the smallest prime strictly less than — — h 2 = 8), q 2 < 31 
(the smallest prime less than 2 2 ' • (9 - 2 + 1) = 32), q 3 < 109, g 4 < 2 8 • 6 = 3 • 2 9 , 
q 5 < 5 • 2 16 , and g 6 < 2 32 • 4 = 2 34 . 

By using the abundancy index function, we can further reduce the bound 
for q 2 . If q 2 > 13, then 

2 = i(N) < /(3 a i3 /3 i7xi9 ,5 23 e 29*iomooo7 r ?iooooooo7 t ) < 

3 • 13 • 17 • 19 • 23 • 29 • 101 • 10007 • 100000007 _ 849216193914429412851 

2 • 12 • 16 • 18 • 22 • 28 • 100 • 10006 • 100000006 ~ 426034693082080051200 ~ 

1.9933029110162608467608441119731 < 2 

so that q 2 < 11. 

We may also try reducing the bound for q 3 . Proceeding in the same manner 
as before, if qi > 3, q 2 > 5 and q 3 > 53, then 

2 = I(N) < 7(3 Q 5 /3 53 >;: 59 <5 61 e 67*101 7 10007 ,7 100000007 t ) < 

3 • 5 • 53 • 59 • 61 • 67 • 101 • 10007 • 100000007 _ 19375328833237423387515 _ 
2 • 4 • 52 • 58 • 60 • 68 • 100 • 10006 • 100000006 ~ 9850131125407832064000 ~ 

1.9670122749188491347643364823597 < 2 

so that g 3 < 47. 

Notice that a major problem with the abundancy index function is that it is 
not capable of determining upper bounds on the prime factors of an OPN beyond 



the smallest three (except in some special cases). This is due to the fact that the 
first three primes could be 3, 5, and 11, in which case the sigma bounds would 
allow an arbitrary number of additional prime divisors, but require no more. 

In a preprint published in the electronic journal INTEGERS in 2003, Nielsen 
improved on Cook's bound by showing that N < 2 4 \ Since qi < N for all i, 
Nielsen's bound is an implicit upper limit on q t . If q t is the special (or Euler) 
prime factor with an exponent of 1 and the other q^s are small, then we can 
say little else about q t . We can, however, give tighter limits for the other prime 
factors. 

Since only one of q t , q t -i can be the special prime, at least one exponent 

1! 

is even. Consequently, qtqt-i 2 < N < 2 4 , so that g t _i < 2 3. Likewise, for 

4* 

1 < i < t, we have qi < 2 2 ( f ~ ^ + 1 . This range can be limited further by 
considering the other prime factors. Hare proved that there are at least 75 total 
primes (not distinct), so we may take the other primes to be as small as possible 
and raise the smallest prime to the appropriate power (and the others to the 2nd 
power). By using such a method, we may be able to reduce the bound by perhaps 
a million, depending on t = ou(N). 

We have thus given explicit formulations (through Griin/Perisastri and Kishore) 
for the upper limits on the smallest 6 prime divisors, which are augmented with 
sigma conditions for the lowest 3. The higher prime divisors are likewise restricted 
(through Nielsen), though not as tightly. 
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Going beyond the results previously mentioned is not easy, considering the 
fact that both Iannucci and Jenkins used proofs based on the divisibility of cy- 
clotomic polynomials F p (x) to find lower bounds for the highest prime factors, 
and that topic is quite hard to follow, to my knowledge. The method used in 
most modern proofs is that of factor /sigma chains. Consider an OPN iV with 
a component 5 4fc+1 . (By component, we mean either a prime power that divides 
N, or simply a factor of N, which may not necessarily be a prime power.) We 
know that er(5) | <r(5 4fc+1 ) for all positive integers k, so we may conclude that 
cr(5) | 2N. But <r(5) = 6, so this is an indication that 3 | N. This illustrates how 
knowledge of a particular prime power leads to knowledge of other prime powers, 
for an OPN, by virtue of the fact that the o function is multiplicative. 

We summarize the results in this subsection (which are taken from Greathouse 
[20J) as follows: 

u(N) = 9 

100000007 < 99 < 2 4 ° 

49 

10007 < qs < 2 3 

49 

101 < <? 7 < 2 5 
23 < 96 < 17179869143 
19 < q 5 < 327673 
13 < (?4 < 1531 
11 < 93 < 47 
5 < q 2 < 11 
3 < qi < 3 



Remark 4.2.2. These bounds, together with the algorithm presented in Subsec- 
tion 4.2.2, can (potentially) prove the conjecture that u(N) > 10, for a general 
OPN N. 
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In the next subsection, we shall discuss some of the author's own results on 
the relationships between the components of an OPN N. 

4.2.4 Relationships Between OPN Components 

Throughout this subsection, we let N = p k m 2 be an OPN with special/Euler 
prime p with p = k = 1 (mod 4) and gcd(p, m) = 1. (Recall from Section 2.4 
that p k is called the Euler's factor of the OPN N.) It will also be useful later 

u>(N) 

to consider the canonical factorization N = J J pi a \ where p 1 < p 2 < ■ ■ ■ < Pt, 

i=i 

t = oj(N) and ai > for all i. 

We begin with some numerical results: 

a(p k ) 5 8 o-(m 2 ) 
Lemma 4.2.3. 1 < < — < — < --^ < 2 

pft 4 5 m z 

Proo/. 



_ p fc +/- i + ...+p 2 + P +i _ i | i | /i\ 2 1 ny- 1 1 (iy >1 , i 

pfc pk p yp J yp J \P J ~ P 



■ , a (P ) 
in particular, — 7— > 1. 



a(p fe ) _ p k+l - 1 



pk pk (j) _ p k (p — 1) p — 1 — - 

But since p k m 2 is an OPN, p is prime and p = 1 (mod 4). This implies that 

1 14 p 15 

p > 5, from which it follows that 1 > 1 = -. Thus, = r < -, 

p 55 p — 1 1 _ £ 4 

P 

o~(v k ^) 5 

and we have — 7— < -• Note that, for k > 1 and prime p, we have: 
p ft 4 
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p + 1 a(p ) p 

< ■ < — — 

p p k p — 1 



Also, since p m is an OPN, 



k\l 



p k 



aim 2 ) 1 
m 2 



2, which implies that 



a(p k ) 2m 2 a(p k ) 2m 2 5 

— r — = —7 — 57. But 1 < — 7— = —7 — ^7 < -. Consequently, we have 
p k aim 1 ) p k aim 1 ) 4 

4 aim 2 ) 1 , , 8 cr(m 2 ) „ ,_, 

5 2m z 5 rrr 

Corollary 4.2.1. If N = p k m 2 is an OPN with Euler's factor p k , then 

p + 1 c(p fc ) P 2 (p — 1) a(m 2 ) ^ 2p 
p — pk p _ 1 p m 2 — p _|_ 

Proof. The proof is similar to that for Lemma 14.2.31 We give here a proof of the 

p 2 (p — 1) 

inequality in the middle. Suppose to the contrary that > . Since 

p—1 p 

p > 5 > 0,p 2 > 2(p — l) 2 . This implies that p 2 — 4p + 2 < 0. This last inequality 
is a contradiction since it implies that p(p — 4) + 2 < 0, whereas p > 5 implies 
that p(p - 4) + 2 > 7. □ 

aim 2 } 

In what follows, we set X = — - — and Y = — . 

pk m 2 

57 aip k ) aim 2 ) 
Lemma 4.2.4. — < + -^r 1 < 3 

20 p k m 2 

Proof. By Lemma EZM l<X<-<-<Y<2. Consider (X - 1)(Y - 1). 

4 5 

This quantity is positive because 1 < X < Y. Thus, 

(X - 1)(Y - 1) = XY - (X + Y) + 1 > 0, 

which implies that X + Y < XY + 1. But XY = 2. Thus, X + Y < 3. Now, 
consider ^X — ^Y — . This quantity is negative because X < ^ < Y. 
Thus, 
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x-^](y--)=xy--(x + y) + — <o, 

A) V V 4 v ; 16 

5 25 
which implies that — (X + Y) > XY + — . But again, XY = 2. Consequently, 

7 (X + Y) > ^, and hence X + Y >^-. □ 
4 lo 20 

Corollary 4 2 2 3 P 2 ~ A P + 2 < a (P k ) + a ( m2 ) < 3p 2 + 2p + 1 
p(p — 1) p k m 2 ~ p(p+l) 

Proof. From Corollary 14.2.11 

p + 1 p 2 (p — 1) 2p 
' < X < — — < '- <Y < 1 



p p — 1 p p + 1 

Consider I X — - ) ( Y — — J . This quantity is nonnegative because 



p J \ p 

p+l 



P 



< X < Y. Thus, 

X-2±T) (V-^Uxy-^(X + Y) + ^>0 
p J \ p J p p 

which implies that - ~*~ (X + Y) < 2 + - ^ = _ Conse- 

p p 2 p 2 

quently, X + Y < ^—t— -^-t- ■ Now, consider ( X — J ( Y — J . This 

p(p + l) V P~ V V P~ V 

quantity is negative because X < — - — < Y. Thus, 

p — 1 

X-^-) ( Y -J P -\=XY-J^{X + Y)+ f P \ 2 <0 

p — i / v p ~ i / p— l (p — l) 2 

p p 2 3p 2 — 4p -|- 2 

which implies that (X + Y) > 2 H — = — — — . Conse- 

p — 1 p 2 — 2p + 1 (p — l) 2 

quently, — - ~[ < X + Y. 
p(p-l) 

Finally, we need to check that, indeed, 

3p 2 - 4p + 2 _ p - 2 p - 1 _ 3p 2 + 2p + 1 

p(p — 1) p(p-l) p(p + l) p(p+l) 

Suppose to the contrary that 



3p 2 - 4p + 2 3p 2 + 2p + 1 
p(p-l) ~~ p(p+l) 

This last inequality implies that 3 ^ > 3 y -, or equivalently, 

p{p - 1) p(p + 1) 

p — 1 p — 2 

> — — -. Since p > 5 > 0, we have (p — l) 2 > (p + l)(p — 2), or 



p(p + 1) p(p - 1) 
equivalently, p 2 — 2p + 1 > p 2 — p — 2, resulting in the contradiction p < 3. 

Hence, we have 

3p 2 -4p + 2 3p 2 + 2p+l 

7 ^ — 7 \ ) 

p(p - 1) p(p+l) 

and we are done. □ 

If we attempt to improve the results of Corollary 14.2.21 using Corollary 14.2. 1[ 
we get the following result: 

Theorem 4.2.1. The series of inequalities 

i(p)< ^ + ^< C7(p) 

p K m z 



with 

Lip) 

and 

U(p) 



3p 2 - 4p + 2 
Pip ~ !) 



3p 2 + 2p + 1 



pip + 1) 

is best possible, for a given Euler prime p = 1 (mod 4) of an OPN N = p k m 2 . 



Proof. From Corollary 14.2. 1} we have: 
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and 



2(P-1) <y< 2p 



p p + 1 

We remark that such bounds for X and Y are best possible by observing that 
k = 1 (mod 4) implies > 1. Adding the left-hand and right-hand inequalities 
give rise to: 

3p-l <x + y< P (3p-1) 



p (p + l)(p-l) 

Comparing this last result with that of Corollary 14.2.21 the result immediately 

follows if we observe that 

f 3p 2 - 4p + 2 3p - 1 1 3p 2 - 4p + 2 

max < - — , > = r — 

{ p(p-l) p J p(p-l) 

and 

. f 3p 2 + 2p + 1 p(3p - 1) 1 _ 3p 2 + 2p + 1 
mm I Kp+T) ' (p+l)(p-l)J ~ p(p + 1) ' 

with both results true when p > 5 (specifically when p is a prime with p = 1 
(mod 4)). □ 

The reader might be tempted to try to improve on the bounds in Lemma 
14.2.41 using Theorem 14.2. 1| but such efforts are rendered futile by the following 
theorem: 



Theorem 4.2.2. The bounds in Lemma 4-2-4 are best possible 



Proof. It suffices to get the minimum value for L(p) and the maximum value for 
U(p) in the interval [5, oo), or if either one cannot be obtained, the greatest lower 
bound for L(p) and the least upper bound for U(p) for the same interval would 
likewise be useful for our purposes here. 
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From basic calculus, we get the first derivatives of L(p), U(p) and determine 
their signs in the interval [5, oo): 

and 

which means that L(p), U{p) are increasing functions of p on the interval [5, oo). 

57 

Hence, L(p) attains its minimum value on that interval at L(5) = — , while U(p) 
has no maximum value on the same interval, but has a least upper bound of 
lim U(p) = 3. 

p^oo 

This confirms our earlier findings that 

57 a{p k ) a(m 2 ) 

on k 1 2~ ' 

20 p K m z 

with the further result that such bounds are best possible. □ 
Remark 4.2.3. Let 

p K a{p K ) 

Using Mathematica, we get the partial derivative: 

d_ _ p^- k (k-kp + p(-l + p k )) (-1 + p k (2p + p k (2 + (-4 + p)p))) 

fyJW )- (_l +p )2(_ 1+p l+fc)2 

which is certainly positive for prime p = 1 (mod 4) and k a fixed positive in- 
teger satisfying k = 1 (mod 4). This means that f(p,k) = X + Y is a strictly 
monotonic increasing function of p, for such primes p and fixed integer k. Lastly, 
lim X + Y = 3. 



Remark 4.2.4. Why did we bother to focus on improving the bounds for X+Y in 
the first place? This is because Joshua Zelinsky, in response to one of the author's 

posts at the Math Forum ( |http://www.mathforum.org/kb/message.jspa?messageID=4140071| fetstart=0), 

said that "[he does not] know if this would be directly useful for proving that no 

[OPNs] exist, [although he is not] in general aware of any sharp bounds [for X+F]. 

Given that there are odd primitive abundant numbers n of the form n = P K M 2 

with P and K congruent to 1 modulo 4 and gcd(P, M) = 1, [he] would be surprised 

if one could substantially improve on these bounds. Any further improvement of 

the lower bound would be equivalent to showing that there are no [OPNs] of the 

form 5m 2 which would be a very major result. Any improvement on the upper 

bound of 3 would have similar implications for all arbitrarily large primes and 

thus [he thinks] would be a very major result. [He's] therefore highly curious as to 

what [the author had] done." In particular, by using Mathematica, if one would 

43 

be able to prove that — < X + Y, then this would imply that p > 5 and we arrive 

15 

at Zelinsky's result that "there are no OPNs of the form 5m 2 " . Likewise, if one 
would be able to derive an upper bound for X + Y smaller than 3, say 2.9995, so 
that X + Y < 2.9995, then this would imply that p < 1999, confirming Zelinsky's 
last assertion. 

Our first hint at one of the relationships between the components of an OPN 
is given by the following result: 

Lemma 4.2.5. If N = p k m 2 is an OPN, then p k ^ m 2 . 

Proof. We will give four (4) proofs of this same result, to illustrate the possible 
approaches to proving similar lemmas: 
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• If p k = m 2 , then necessarily u(p k ) = u(m 2 ). But u(p k ) = 1 < 8 < u(m 2 ), 
where the last inequality is due to Nielsen. 

• Suppose p k = m 2 . This can be rewritten as p ■ p k ~ l = m 2 , which implies 
that p | m 2 since p is a prime. This contradicts gcd(p, m) = 1. 

• Assume p k = m 2 . Then N = p 2k is an OPN. This contradicts the fact that 
prime powers are deficient. 

• Let p k = m 2 . As before, N = p 2k is an OPN. This implies that <j(N) = 
a(p 2k ) = l+p + p 2 + ...+p 2k - 1 +p 2k = (2k + 1) (mod 4) = 3 (mod 4) (since 
p = k = 1 (mod 4)). But, since N is an OPN, a(N) = 2N. The parity of 
LHS and RHS of the equation do not match, a contradiction. 

□ 

By Lemma I4.2.5[ either p k < m 2 or p k > m 2 . 

We now assign the values of the following fractions to the indicated variables, 
for ease of use later on: 

g(p fc ) 
P 1 = — ~ 

a(p k ) 

P2 = — 

m z 

cr(m 2 ) 

Pi = — 

m z 

a(m 2 ) 

pK 

5 g 

From Lemma 14.2.31 we have 1 < p\ < — < — < H\ < 2. Note that pifii = p 2/ u 2 = 

4 5 

2. Also, from Lemma I4.2.5[ we get p\ ^ p 2 and /ii ^ /i 2 - 
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The following lemma is the basis for the assertion that "Squares cannot be 
perfect" , and will be extremely useful here: 

Lemma 4.2.6. Let A be a positive integer. Then c(A 2 ) is odd. 

R 

Proof. Let A = J^Jft^ be the canonical factorization of A, where R = u(A). 
3=1 

Then 

R Ft 

-(^ 2 )=-(n^)=n^ 2ft ), 

3=1 3=1 

since a is multiplicative. But 

fi^) = n i 1 + * + ^ 2 + • • • + ^ 2ft_i + 

3=1 3 = 1 

But this last product is odd regardless of whether the q^s are odd or even, i. e. re- 
gardless of whether A is odd or even. Consequently, o~(A 2 ) is odd. □ 

A reasoning similar to the proof for Lemma 14.2.61 gives us the following 
lemma: 

Lemma 4.2.7. Let N = p k m 2 be an OPN with Euler's factor p h . Then 

a(p k ) ^ a(m 2 ). 

Proof. By Lemma |4.2.6[ a{m 2 ) is odd. If we could show that a{p k ) is even, then 
we are done. To this end, notice that 

a(p k ) = 1 + p + p 2 + ... + p k = (k + 1) (mod 4) 

since p = 1 (mod 4). But k = 1 (mod 4). This means that a(p k ) = 2 (mod 4), 
i. e. a(p k ) is divisible by 2 but not by 4. □ 
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a(j) k ^) aim* 2 } 

From Lemma 14.2.71 we get at once the following: p 1 = — — ^ — = p 2 

pk pk 

. a(p k ) aim 2 ) ... 

and p 2 = 5— f 1 7T- = p\- Also, by a simple parity comparison, we get: 

m 2 m 1 

a(p k ) . aim 2 ) 
p 2 = 7T- t 1— = A*2- Lastly, p 2 t 1 and p 2 t 2 - 

rrr p K 

From the equation a(p k )a(m 2 ) = 2p k m 2 and Example 14.1.41 we know that 
p k I a(m 2 ). This means that p 2 G Z + . Suppose that p 2 = 1. This implies that 
p 2 = 2, and therefore, a(m 2 ) = p k and a(p k ) = 2m 2 . However, according to the 
paper titled "Some New Results on Odd Perfect Numbers" by G. G. Dandapat, 
J. L. Hunsucker and Carl Pomerance: No OPN satisfies a(p k ) = 2m 2 , a{m 2 ) = p k . 
[T5] This result implies that p 2 ^ 2 and p 2 ^ 1. But p 2 7^ 1 and p 2 7^ 2. Since 
p 2 G Z + , we then have p 2 > 3. (Note that p 2 must be odd.) Consequently, we 
have the series of inequalities: 

2 5 8 
0<p 2 <^<l<Pi<T<3</ii<2<3</i 2 . 

3 4 5 

In particular, we get the inequalities p k < a(p k ) < -m 2 and ^ ] = ^ 1 ^ 2 < 

3 a{m z ) pi + p 2 

5 

12' 

57 

Recall that, from Lemma 14.2.41 — < pi + pi < 3. Consider (p 2 — 3)(p 2 — 3). 

20 

This quantity is nonpositive because p 2 < 3 < p 2 . But (p 2 — 3)(/i 2 — 3) < implies 



that p 2 p 2 — 3(p 2 + ^2) + 9 < 0, which means that 11 < 3(p 2 + P2) since p 2 p 2 = 2. 
Consequently, we have the series 
In particular, p\ + p\ 7^ p 2 + /i 2 . 



57 11 
Consequently, we have the series of inequalities — < p\ + p\ < 3 < — < p 2 + p 2 . 

20 3 



We summarize our results from the preceding paragraphs in the theorem 
that follows: 
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Theorem 4.2.3. 

2 5 8 
0<p 2 <-<l<Pi<7<-<Pi<2<3</x 2 

3 4 5 

and 

57 11 

— < pi + in < 3 < — < p 2 + H2 

Remark 4.2.5. We remark that the results here were motivated by the initial 
finding that I(p k ) < I(m 2 ), i. e. p\ < p,i. We prove here too that I(p k ) < 
I(m). We start with: For all positive integers a and b, a(ab) < a(a)a(b) with 
equality occurring if and only if gcd(a, b) = 1. (For a proof, we refer the interested 
reader to standard graduate textbooks in number theory.) It is evident from this 
statement that for any positive integer x > 1, I(x 2 ) < (I(x)) 2 . In particular, from 



8 2 v 10 

Lemma 14.2.31 we have - < I{m 2 ) < (/(m)) 2 , which implies that < I(m). 

5 ' 5 



But I(p k ) < 1.25 (again from Lemma 12X31 . and ^— rj 1.26491106406735. 

5 

Consequently, I{p k ) < I{m). (Note that gcd(p, m) = 1.) This should motivate 

2 

the succeeding discussion, which attempts to improve the result p k < -m 2 to 
p k < m, where again iV = p k m 2 is an OPN with Euler's factor p k . 



We now attempt to obtain the improvement mentioned in Remark 14.2.51 
Since proper factors of a perfect number are deficient, we have I(m) < 2. Conse- 



quently, we have the bounds < I(m) < 2. Likewise, since p k m is a proper 

5 

factor of iV which is perfect, and from Lemma [4.2.31 we have 1 < I(p k ), hence the 
following must be true: 

Lemma 4.2.8. Let N = p k m 2 be an OPN with Euler's factor p k . Then 



5 m p k 
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Now, as before, let 



a(p k ) 

pie 

a(p k ) 

ps = 

m 

aim) 

A*3 — 

m 

a(m) 

= — 7— 



,i 2V10 , 2V10 

From the preceding results, we get 1 H < pi + p 3 < 3, where 1 H ~ 

5 5 

2.26491 and the rightmost inequality is obtained via a method similar to that 

used for Lemma 14.2.41 Also, from Lemma 14.2.81 and Lemma 13.2.11 (Arithmetic 

2 v^lOOO 

Mean- Geometric Mean Inequality) , we get the lower bound < P3 + A*4, 

5 

where « 2.24937. 



We observe that P3 7^ 1 since <r(p fc ) = 2 (mod 4) while m = 1 (mod 2) since 
TV is an OPN. Hence, we need to consider two (2) separate cases: 

Case 1: p 3 < 1. Multiplying both sides of this inequality by p 4 > 0, we get 



P3/U4 < fi4. But from Lemma I4.2.8[ we have — - — < p 3 p 4 . Therefore: p 3 < 1 < 



2\/l0 a(v k ) aim) 

1.26491106406735 « — — < p 4 . (In particular, p 3 = ^ = ^ 4 .) Note 

5 m p k 



u 
p*- 

that, from Lemma [4.2.31 — < p 3 < 1, which implies that p h < m. 

m 



Case 2: 1 < p 3 . Similar to what we did in Case 1, we get from Lemma 

, , , , (h/A)p k 

14.2.81 u 4 < 2. Note that, from Lemma 1472731 1 < p 3 < ±-!—L — which implies 

m 

4 

that p k > —m. 
y 5 
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We now claim that p 3 7^ p 4 also holds in Case 2. For suppose to the 
contrary that 1 < p 3 = p 4 < 2. Since these two ratios are rational numbers 
between two consecutive integers, this implies that m \ a(p k ) and p k \ a(m). But 
P3 = P4 =>■ p k <y(p k ) = ma(m), which, together with gcd(p, m) = gcd(p fe ,m) = 1, 
results to a contradiction. This proves our claim. Hence, we need to consider two 
(2) further subcases: 

Subcase 2.1: 1 < p 3 < p 4 < 2 



mcr(m) 



^^^^ P^r 

m p K \p K my {p K m) z 



1 (j(p fe ) ^ 1 cr(m) 1 1 a(p k ) ^ 1 a(m) ^ 2 
m m 2 m 2 p k 
p 2k < 2m 2 p k < V2m 



m 2 pk p2k m m 2 m 2 pk p2k m p2k 



Thus, for Subcase 2.1, we have: —m < p k < \p2m. (Note that it may still 

5 

be possible to prove either p k < m or m < p k in this subcase. We just need to 
further develop the methods used.) 

Note that we can improve the bounds for p 3 and p 4 in this subcase to: 

r- , ^lOOO 
1 < P3 < V 2 and — - — < /i 4 < 2, 

5 

and that we can derive the upper bound 

p 3 + P4 < 3. 

Subcase 2.2: p 4 < p 3 , 1 < p 3 and p 4 < 2 

cr(m) <r(p fc ) k / k\ maim) p k a(p k ) 

pK m \p K m,y \p K my 



Ill 



1 a{m) 1 a(p k ) 1 a(m) 



p2k m m 2 pk m 2 m 
-zr < — =>- vri < p lK =4> m < p k 

It is here that the author's original conjecture that p k < m in all cases is 
disproved. 

Note also the following improvements to the bounds for p 3 and /i 4 in this 

subcase: — - — < p 3 and p 4 < V A. 
5 

At this point, the author would like to set the following goals to treat this 
Subcase 2.2 further: 

• Obtain an upper bound for p 3 . 

• Obtain a lower bound for p 4 . 

• Obtain an upper bound for p 3 + p 4 . 

We summarize our results in the following theorem: 

Theorem 4.2.4. Let N = p k m 2 be an OPN with Euler's factor p k . Then p 3 ^ p A , 
and the following statements hold: 

• If p 3 < 1, then p k < m. 

• Suppose that 1 < p 3 . 

- If p 3 < Pa, then \m < p k < \plm. 

5 

— If pi < p 3; then m < p k . 
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2 

We now state and prove here the generalization to P2 < - mentioned in 

3 

Section 3.2: 

Theorem 4.2.5. Let N = p^ be the canonical factorization of an OPN N, 

t=i 

where Pi < P2 < ■ ■ ■ < Pt (ire primes, t = u(N) and a* > for all i. Then 
2 N 

o-{pi ai ) <-— for all i. 

3pi a% 

Proof. Let A" = pi ai M for a particular i. Since Pi ai \\N and A" is an OPN, then 
a(pi ai )a(M) = 2pi ai M. From Example 14.1.41 we know that pi ai \ a(M) and we 
have a(M) = hpi ai for some positive integer h. Assume h — 1. Then a(M) = Pi a \ 
forcing a(pi ai ) = 2M. Since A^ is an OPN, p^ is odd, whereupon we have an odd 
cx,i by considering parity conditions from the last equation. But this means that 
Pi ai is the Euler's factor of N, and we have Pi ai = p k and M = m 2 . Consequently, 
a(m 2 ) = a(M) = pi ai = p k , which contradicts the fact that /i 2 > 3. Now suppose 
that h = 2. Then we have the equations a(M) = 2pi ai and o~{pi ai ) = M. (Note 
that, since M is odd, «j must be even.) Applying the a function to both sides of 
the last equation, we get a(a{pi ai )) = a(M) = 2pi ai , which means that pi ai is an 
odd superperfect number. But Kanold [66J showed that odd superperfect numbers 
must be perfect squares (no contradiction at this point, since a« is even), and 
Suryanarayana [67] showed in 1973 that "There is no odd super perfect number 

of the form p 2an (where p is prime). Thus h = — — - > 3, whereupon we have 

Pi ai 

2 2 N 

the result a(pi ai ) < -M = for the chosen i. Since % was arbitrary, we have 

3 3 Pi a * 

proved our claim in this theorem. □ 



The following corollary is a direct consequence of Theorem 14.2.51 



Corollary 4.2.3. Let N be an OPN with r = u(N) distinct prime factors. Then 




In the next section, we attempt to "count" the number of OPNs by trying 
to establish a bijective map between OPNs and points on a certain hyperbolic 
arc. We prove there that such a mapping (which is based on the concept of the 
abundancy index) is neither surjective nor injective, utilizing results on solitary 
numbers in the process. 



4.3 "Counting" the Number of OPNs 

In this section, we will be disproving the following conjecture: 

Conjecture 4.3.1. For each N = p k m 2 an OPN with N > 10 300 , there corre- 
sponds exactly one ordered pair of rational numbers \ ^ , ° m ] lying in the 

\ p k m 2 J 

a(p k ) 5 8 aim 2 ) , 57 a(p k ) aim 2 ) 

region 1 < < - - < --^ < 2, and — < + --^ < 3, and 

pK 4 5 m z 20 p K m z 

vice-versa. 

a(p k ) 

We begin our disproof by observing that, for prime p, — r— is a decreasing 



function of p (for constant k) and is also an increasing function of k (for constant 

Pl kl P2 k 

(1) p\ 7^ P2, k\ = k 2 and (2) p\ = p 2 , k\ ^ k 2 . To show that the same inequality 



p). These observations imply that \ — ^ \ — for the following cases: 



holds for the case (3) p\ ^ p 2 ,k\ ^ k 2 , we proceed as follows: Suppose to the 

Pl kl P 2 k 



contrary that : — = : — but p\ ^ p 2 ,k\ ^ k 2 . Then 



p 2 k2 {p 2 - l)( Pl kl+1 - 1) = Pl k ^ Pl - l){p 2 k2+1 - 1). 
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Since pi and p 2 are distinct (odd) primes, gcd(pi,p 2 ) = 1 which implies that 

Pi kl | (P2 - l)(pi hl+1 - 1) and p 2 k2 | (pi - \){p 2 k ' 2+l - 1). 
Now, let us compute gcd(p k ,p k+1 — 1) using the Euclidean Algorithm: 

pk+l _ i _ p . p k _ ]_ _ ^ _ 1) . pk _|_ (pfe _ 1) 

p fc = 1 ■ (p fe — 1) + 1 — > last nonzero remainder 

pfc - 1 = (p* - 1) . 1 + 

Consequently, gcd(p fc ,p fe+1 — 1) = 1. Hence, we have that 

Pi kl I (P2 — 1) and p 2 k2 | (pi — 1). 

Since p\ and p 2 are primes, this means that p\p 2 \ (pi — l)(p 2 _ 1), which further 
implies that p\p 2 < (pi — l)(p 2 — 1) = pip 2 — (pi + p 2 ) + 1, resulting in the 
contradiction p 1 + p 2 < 1. 

Remark 4.3.1. An alternative proof of the results in the preceding paragraph 
may be obtained by using the fact that prime powers are solitary (i. e. I{p k ) = 
I(x) has the sole solution x = p k ). 

<r(p k ) (j(Tn?) 

Let X = — 7— and Y = — . It is straightforward to observe that, since 

pk m 2 

the abundancy index is an arithmetic function, then for each iV = p k m 2 an OPN 
(with iV > 10 300 ), there corresponds exactly one ordered pair of rational numbers 
(X,Y) lying in the hyperbolic arc XY = 2 bounded as follows: 1 < X < 1.25, 
1.6 < Y < 2, and 2.85 < X + Y < 3. (Note that these bounds are the same ones 
obtained in Subsection 4.2.4.) 
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We now disprove the backward direction of Conjecture 14.3.11 We do this by 

} (tItTI 2 ^] 

showing that the mapping X = — - — and Y = — is neither surjective nor 

injective in the specified region. 



(X, Y) is not surjective. We prove this claim by producing a rational 

a(p k ) 

point (X ,Y ) lying in the specified region, and which satisfies X ^ — — for 

pk 

all primes p and positive integers k. It suffices to consider the example Xq = 

where p and q are primes satisfying 5 < p < q. Notice that 1 < Xq = 

pq 

= (l + -)(l + -)<--< 1-2468 < 1.25. Now, by setting 
pq v P/ \ Qj 711 

2 

Y = — , the other two inequalities for Y and X + Y would follow. Thus, we 

now have a rational point (Xo,Yo) in the specified region, and which, by a brief 

a(p k ) 

inspection, satisfies X — 7— for all primes p and positive integers k (since 

pk 

prime powers are solitary). Consequently, the mapping defined in the backward 
direction of Conjecture 14.3. II is not surjective. 



Remark 4.3.2. Since the mapping is not onto, there are rational points in the 
specified region which do not correspond to any OPN. 

(X, Y) is not injective. It suffices to construct two distinct OPNs 
^1 = Pi kl mi 2 and iV 2 = P2 k2 m2 2 (assuming at least two such numbers exist) that 
correspond to the same rational point (X, Y). Since it cannot be the case that 
Pi kl 7^ p2 k2 ,mi 2 = m2 2 , we consider the scenario pi kl = p2 k2 ,rrii 2 7^ ni2 2 . Thus, we 
want to produce a pair (mi, 777,2) satisfying I(mi 2 ) = J(m 2 2 ). (A computer check 
by a foreign professor using Maple produced no examples for this equation in the 
range 1 < rrii < m 2 < 300000. But then again, in pure mathematics, absence 



of evidence is not evidence of absence.) Now, from the inequalities p k < m? and 

N = p h m 2 > 10 300 , we get m 2 > 10 150 . A nonconstructive approach to finding 

a solution to I(mi 2 ) = J(m 2 2 ) would then be to consider 10 150 < m{ 2 < m 2 2 

and Erdos' result that "The number of solutions of 1(a) = 1(b) satisfying 

g 

a < b < x equals Cx + o(x) where C > — ." ([IT], [2]) (Note that C here 

is the same as the (natural) density of friendly integers.) Given Erdos' result 

g 

then, this means that eventually, as mi — > oo, there will be at least ^ 15 °Y^f 
solutions (1711,1712) to /(mi 2 ) = I(m2 2 ), a number which is obviously greater than 
1. This finding, though nonconstructive, still proves that the mapping defined in 
the backward direction of Conjecture 14.3.11 is not injective. 

Thus, we have failed our attempt to "count" the number of OPNs by showing 
that our conjectured correspondence is not actually a bijection. However, this 
should not prevent future researchers from conceptualizing other correspondences 
that may potentially be bijections. The author is nonetheless convinced that, 
ultimately, such bijections would have to make use of the concept of the abundancy 
index one way or another. 



Chapter 5 



Analysis and Recommendations 



Up to this point, the OPN problem is still open. We have discussed some of the 

old as well as new approaches to solving this problem, with the factor chain 

approach taking center stage in mainstream OPN research. But it appears that, 

for the ultimate solution of the problem using such approaches, we need fresh 

ideas and what may be called a "paradigm shift". For example, the index/outlaw 

status of the fraction — ^- — - , which may potentially disprove the OPN Conjec- 

V 

ture, might probably require concepts from other disparate but related fields of 
mathematics, such as real analysis. But then again, this is just an educated guess. 

Two open problems necessitating further investigations would come to mind, 
if the reader had read and fully understood the contents of this thesis: 

• Prove or disprove: If N — p k m 2 is an OPN with Euler's factor p k , then 
p k < m. 

• There exists a bijection (X, Y) from the set of OPNs to the hyperbolic arc 
XY = 2 lying in the region 1 < X < 2, 1 < Y < 2, 2^/2 < X + Y < 3, 
and this bijection makes use of the concept of abundancy index to define 



the mappings X and Y. 

The underlying motivation for these open problems have been described in suffi- 
cient detail in Chapter 4 and should set the mood for extensive investigation into 
their corresponding solutions. 

Lastly, we warn future researchers who would be interested in pursuing this 
topic that while a multitude of evidence certainly suggests that no OPNs exist, 
neither heuristic nor density arguments but only a (mathematical) proof of 
their nonexistence could lay the problem to a "conclusive" rest. The author has 
tried, and although he failed like the others, eventually somebody will give a final 
end to this problem that had defied solution for more than three centuries. 
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